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Abstract. In this paper we present a new unified theoretical framework that 
describes the full dynamics of quantum computation. Our formulation allows 
any questions pertaining to the physical behavior of a quantum computer to 
be framed, and in principle, answered. We refer to the central organizing prin- 
ciple developed in this paper, on which our theoretical structure is based, as 
the Quantum Computer Condition (QCC), a rigorous mathematical statement 
that connects the irreversible dynamics of the quantum computing machine, 
with the reversible operations that comprise the quantum computation in- 
tended to be carried out by the quantum computing machine. Armed with the 
QCC, we derive a powerful result that we call the Encoding No-Go Theorem. 
This theorem gives a precise mathematical statement of the conditions under 
which fault-tolerant quantum computation becomes impossible in the pres- 
ence of dissipation and/or decohcrcncc. In connection with this theorem, we 
explicitly calculate a universal critical damping value for fault-tolerant quan- 
tum computation. In addition we show that the recently-discovered approach 
to quantum error correction known as "operator quantum error-correction" is 
a special case of our more general formulation. Our approach furnishes what 
we will refer to as "operator quantum fault-tolerance." In particular, we show 
how the QCC allows one to derive error thresholds for fault tolerance in a 
completely general context. We prove the existence of solutions to a class of 
time-dependent generalizations of the Lindblad equation. Using the QCC, we 
also show that the seemingly different circuit, graph- (including cluster-) state, 
and adiabatic paradigms for quantum computing are in fact all manifestations 
of a single, universal paradigm for all physical quantum computation. 
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1. Introduction 

The promise inherent in quantum computing has stimulated a tremendous ex- 
plosion of interest in the research community. Voluminous research has been carried 
out directed to many different problems associated with the development and prop- 
erties of quantum computational algorithms. Parallel to these efforts, substantial 
investigations have been devoted to the problems associated with the development 
of actual quantum computing machines. A rigorous and fully general theory that 
connects quantum computing algorithms and quantum computing machines would 
be of considerable value. 

In this paper we present a new unified theoretical framework that describes 
the full dynamics of quantum computation. Our formulation allows any questions 
pertaining to the physical behavior of a quantum computer to be framed, and 
in principle, answered. We refer to the central organizing principle developed in 
this paper, on which our theoretical structure is based, as the Quantum Computer 
Condition (QCC), a rigorous mathematical statement that connects the irreversible 
dynamics of the quantum computing machine, with the reversible operations that 
comprise the quantum computation intended to be carried out by the quantum 
computing machine. 

The actual dynamics of the system that we intend to use as a practical quantum 
computing machine are those of an open quantum mechanical system, burdened 
with various dissipative and/or decoherence effects. The QCC provides a set of 
mathematical constraints that must be satisfied by a physical system if we intend 
to use that system as a quantum computing machine. 

Armed with the QCC, we derive a powerful result that we call the Encoding 
No-Go Theorem. The Encoding No-Go theorem gives a precise mathematical state- 
ment of the conditions under which fault-tolerant quantum computation becomes 
impossible in the presence of dissipation and/or decoherence. We provide a rigorous 
definition of damping, which includes the phenomena of dissipation and decoher- 
ence, and explicitly calculate a universal critical damping value for fault-tolerant 
quantum computation. This fundamental theorem has deep formal significance. 
Moreover, it also furnishes criteria for solving diverse problems associated to ac- 
tual physical quantum computer realizations, such as determining which practical 
design choices for quantum computing machines are not viable. 
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In addition we show that the recently-discovered approach to quantum error 
correction known as "operator quantum error-correction" (OQEC) is actually a 
special case of our general formulation. Our approach furnishes what we will refer 
to as "operator quantum fault-tolerance" (OQFT). In particular, we show how the 
QCC allows one to derive error thresholds for fault tolerance in a completely general 
context. 

In this paper we define the concept of a quantum component, which allows us 
to study realistic implementations of quantum computers, in which decoherence 
and/or dissipative effects are present, using a dynamical equation of motion suit- 
able for describing an open quantum system. By using the QCC, we are able to 
reconcile the apparent contradiction between: (1) the fact that quantum compu- 
tations are specified by unitary transformations, the associated dynamics of which 
are intrinsically reversible, and (2) the fact that quantum computers, qua practical 
machines, are inevitably characterized by irreversible dynamics. The reconcilia- 
tion suggests an analogy with the fluctuation-dissipation theorem, which relates 
irreversible dynamics to equilibrium properties in a large class of physical systems. 

In this paper we present an existence proof for fundamental solutions to useful 
classes of time- dependent generalizations of the Lindblad equation. This provides 
a useful tool in analyzing a wide variety of open quantum mechanical systems. 

Our framework is sufficiently general to encompass, and describe in a unified 
manner, the currently-known "paradigms" for quantum computation, including 
the circuit-based ("two-way computing") paradigm, the graph state-based ("one-way 
computing") paradigm and the adiabatic quantum computer paradigm. Using the 
QCC, we show by explicit construction that these seemingly different paradigms are 
in fact all manifestations of a single, universal paradigm for all physical quantum 
computation. 1 

2. The Quantum Computer Condition 

2.1. Introduction. In this section we present the Quantum Computer Con- 
dition, a rigorous mathematical statement of the constraints that determine the 
viability of any practical quantum computing machine. To achieve the goal of 
practical quantum computation we must produce an actual physical device that 
implements a predetermined unitary operator U acting on some Hilbert space. The 
Quantum Computer Condition relates the unitary operator representing a quantum 
computation to the actual physical device intended to perform that computation. 

The specification of U defines ideally the quantum computation to be performed 
by the quantum computing machine. Gencrically, the result of a quantum compu- 
tation, U, is then used to carry out the probabilistic evaluation of some classical 
function. The complete quantum computation comprises a number of elements, 
including 

• Preparation of a quantum state for initialization. 

• Measurement of a quantum state for readout. 



In the particular case of the graph state-based paradigm (which includes cluster state-based 
models), we not only show that the paradigm is a manifestation of the unifying picture provided by 
the QCC, but also introduce a definition of graph state-based quantum computers that generalizes 
the graph state models previously defined in the literature. 
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• Various tasks that can be performed by classical computers, such as pre- 
processing of the data, or postprocessing of the output into some humanly 
comprehensible form. 

However, the above list of elements are not what is "important" about quantum 
computers. Rather: 

• The distinctive element of quantum computation is the "ability to perform 
quantum gates" (c.f. [23], §4.6). 

Mathematically, a quantum gate is a unitary (hence reversible) operator U 
acting on a Hubert space. The formally defined "gates," as such, are not "devices." 
They arc concepts: they don't implement themselves. A machine is required to 
physically implement the abstractly defined unitary transformation. The actual, 
physical computing device intended to implement the transformation is described 
mathematically by a completely positive trace-preserving map, P, that transforms 
the input state to the output state. We will refer to a physically realizable device 
intended to implement an ideal quantum computation as a quantum component. In 
this paper we study realistic quantum components, in which decoherence and/or 
dissipative effects are present, using a dynamical equation of motion suitable for 
describing an open quantum system. 2 We must reconcile the fact, and apparent 
paradox, that a non-reversible mapping, P. is used to "implement" a reversible one, 
U. 

2.2. The Motivation of the Quantum Computer Condition. Mathe- 
matically, a quantum computation is a unitary operator U in the unitary group 
of a Hilbert space. A quantum component is described by a completely positive 
trace-preserving map P which maps the set of trace class operators on the Hilbert 
space to itself. The map P accounts for decoherence and dissipation, as well as 
unitary evolution. We will subsequently discuss in more detail the actual form for 
P. In our analysis we will consider the action of P on density matrices p i— > P ■ p 
rather than on state vectors (and correspondingly the action of U on density matri- 
ces p i — ► UpW, rather than the action of U on state vectors). This is because, due 
to the presence of decoherence and/or dissipation, our system will almost always 
evolve into a mixed state, which can only be described by a density matrix p. 3 

In order to motivate the Quantum Computer Condition (QCC), let us first con- 
sider the abstractly-defined quantum computation itself, prescribed by the unitary 



In order to analyze the effects of dissipation and/or decoherence one must use some method 
of approximating the dynamics of the degrees-of-freedom comprising the rest of the universe 
"outside of the quantum computer. This is of course because the complete, detailed, exact 
analytical solution to the Schrodinger equation of the universe, for all degrees-of-freedom, is not 
known. One reasonable approach is to construct a Lindblad-type equation, based on a presumption 
of underlying Markovian dynamics, in which environment degrees-of-freedom are traced over in 
such a way as to result in a first-order (in time) differential equation. In this paper, for definiteness, 
we utilize a generalized Lindblad-type equation to describe the environment: this is used merely 
in order to exemplify how one may take into account the effects of dissipation and/or decoherence. 
However, most of the results in our paper, including the crucial Encoding No-Go Theorem, are 
independent of this choice, and in particular are independent of the assumption of underlying 
Markovian dynamics. 

''Other reasons for utilizing density matrices rather than state vectors include the generic 
importance in quantum information theory of trace-preserving completely positive maps (which 
restrict to transformations on density matrices), and the useful algebraic and analytic properties 
of density matrices (density matrices for instance form a weak-* compact convex set). 
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operator U. This is assumed to be given, and is represented by 

(1) abstract computation: UpW . 

The action of the quantum component intended to effect the computation is repre- 
sented by 

(2) practical implementation: P ■ p . 

Motivated by the notion of having the machine implement the computation, if we 
were to require that the identity 

(3) P ■ p = UpU^ 

hold for all density states p, then the action of P would in fact be identical to the 
action of the unitary operator. This would imply that P preserves von Neumann 
entropy ( |32j . §5.3), and hence actually models a system with neither dccoherence 
nor dissipation, which is not the case for a practical quantum computing machine. 
Thus, equation Q cannot furnish the correct constraints for realistic quantum 
computation. We will accordingly refer to equation J2Jl as the ersatz quantum 
computer condition (£QCC). 4 

More realistically, taking into account the inevitable presence of decoherence, 
we can require that J3J hold for some restricted set of density states. In this case, 
the solution set will correspond to decoherence-free subspaccs. In order to analyze 
this, we must carefully distinguish between the two different Hilbert spaces that 
arise in this problem. The abstract quantum computation is defined on the Hilbert 
space of logical quantum states, -ffiogicai, so that we have 

(4) U : i?logical — * ^logical- 

In contrast, the presence of decoherence (which affects the actual device) necessi- 
tates that the completely positive trace-preserving map P (which represents the 
actual device) is associated to a different Hilbert space, ff comp , the states of which 
are referred to as computational quantum states. The decoherence-free subspace is 
contained within H comp . (The specific decoherence is accounted for in the explicit 
form of P). As noted above, a consequence of the decoherence is that P operates 
on density matrices rather than on state vectors. Letting T(_ff comp ) be the Banach 
space of trace-class operators on H comp , we have 

(5) P : T(H comp ) -> T(F comp ). 

In order to replace £QCC (equation Q) with an equation that properly incor- 
porates decoherence effects so that it can be used to determine decoherence-free 
subspaces, we must introduce suitable encoding and decoding maps that connect 
the relevant Hilbert spaces. We can hope to find an encoding operator M C nc de- 
fined on a space of logical inputs and a decoding operator Afdcc with values in a 
space of logical outputs, the actions of which are given by (here T(ifi og icai) is the 
Banach space of trace-class operators on -ffi og icai, analogous to T(7J comp )) 

(6) Menc ■ T(ff logical) — > T(i? comp ) 



Although the £QCC does not describe practical quantum computing machines, we note that 
it can be shown that in the finite-dimensional case, the set of p which satisfy J3J is an algebra 
2lp u which depends on both P and U . The details of how one explicitly obtains the algebra 2lp y 
of solutions to J2J are given in Appendix Icl 
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and 




A^dcc : T(JJ C omp) — * T(-fflogical), 



for all logical inputs p. We will refer to equation ((SJ as the "encoded ersatz quan- 
tum computer condition (efQCC). 5 The existence of the encoding and decoding 
maps Alone and .Mdcc is a consequence of the presumed existence of an associ- 
ated decoherence-free subspace of H comp , of dimension greater than or equal to the 
dimension of H logical- 

The meaning of e£ QCC given in eq.JHJ) is as follows. Given a chosen quantum 
computation, U, we wish to construct a physical "machine," P, that implements U. 
In order to do this we must find encoding and decoding maps Ai cnc and A^dec such 
that the equation is satisfied for all p. This is a crucial difference between eqs.© 
and @ : requiring that eq. © holds for all p implies a machine that preserves von 
Neumann entropy, does not dissipate heat and does not decohere, and thus does 
not describe a practical quantum computing device. In contrast, eq.(jHJ holds for all 
states p, but does so by making use of the encoding and decoding operators to map 
to a decoherence-free subspace. The e£QCC (eq.ljHJ) formally holds for all density 
states p. similar to eq.lO, but the use of the encoding and decoding maps in eq.JHJ) 
effectively confines the solution to a restricted set in H comp . 

However, eq.(JHJ) docs not in general provide an acceptable condition to connect 
the dynamics of a practical quantum computing device to the constraints implied 
by the unitary operator U that defines the abstract quantum computation. This 
is because the formulation presented by © does not address situations in which 
residual errors cannot be completely eliminated, even with the use of decoherence- 
free subspaces, and/or other error correction methods |20| . |28| . 

To recapitulate, eq.J2J) describes a quantum computer that performs the re- 
quired computation U, but only if the implementing device, described by the com- 
pletely positive map P, neither dissipates heat nor decoheres. We thus reject this 
expression as a viable quantum computer condition because it describes a system 
that is effectively impossible to achieve. In contrast, cq.(|SJ) describes a quantum 
computer that performs the required computation U, but only if the device de- 
scribed by the completely positive map P implements the required decoherence-free 
subspace with absolutely no residual errors. This does not provide a sufficiently 
general formulation: we need to consider situations in which residual errors cannot 
be completely eliminated. 

2.3. The Presentation of the Quantum Computer Condition. We wish 
to allow for the likelihood that, even if a decoherence-free subspace can be found, 
and even if error correction procedures are applied, there will still be residual errors 
characterizing the operation of the quantum computer. Such a situation may arise 
even if error correction is correctly applied: for instance, in the application of 
concatenated error codes, one iterates the concatenation process until the error 
probability is reduced to a value that is deemed "acceptable" |28j , |24j . This final 
error probability, though small, is not exactly zero. The important point is that 

^Note that no encoding map is required on the right-hand side of equation JSJ since the 
unitary U by definition acts on -Hi og i ca i. 
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it is prudent to write down our quantum computer condition so as to reflect the 
inevitable survival of some amount of residual error. 

In order to quantify the degree to which the actual computational device, repre- 
sented by P, cannot exactly (because of residual error) implement the ideal quantum 
computation, represented by U, we consider the following difference (c/cq.©) 

(9) M dec (P ■ (M enc (p))) - UpUl 

We now compute for this difference a suitable norm on matrices (this norm is made 
more precise below), as 

(10) \\M dcc (P ■ (M cn c(p))) - UpU^\\. 

This quantity is of fundamental importance: it is a measure of the inaccuracy of 
the implementation of U by P. It tells us how well a practical quantum computing 
device actually implements an ideally-defined quantum computation. We will refer 
to the scalar quantity given by l(l()|) as the implementation inaccuracy associated to 
the pair U and P. In connection with this, we introduce a parameter, a, to specify 
the maximum tolerable implementation inaccuracy, so that we have 

(11) \\M dec (P ■ (M cnc (p))) - UpU*\\ < a. 

2.3.1. Formal Statement of the Quantum Computer Condition (QCC). Moti- 
vated by the above considerations, we now introduce and formally define an inequal- 
ity of fundamental importance in the theory of quantum computation that we will 
refer to as the Quantum Computer Condition. We will formally designate this con- 
dition by QCC(P, U, Alone, M. dec 7 ex). In the following, we will impose no constraint 
on the dimensionality of the Hilbert space, and in particular we allow Hilbert spaces 
of infinite dimensions. Let U be a unitary on -ffi og icai and P be a trace-preserving 
completely positive map on T(iJ comp ). Let M cnc : T(7Ji ogical ) — > T(H comp ) and 
A^dec : T(7J comp ) — > T(_ffi og j ca i) be completely-positive, trace-preserving encoding 
and decoding maps (with no further restrictions of any kind on the encoding and 
decoding maps). The Quantum Computer Condition, QCC(P, U, M. enc , Mdec, «), 
holds iff for all density matrices p £ T(Zfi og i ca i), we have 

(12) \\M dcc (P • (AW?))) - UpU^h < a, 

where T(i? comp ) and T(iJi ogica i) are the Banach spaces of trace-class operators 
on -ffcomp and i?i og icai, respectively, and || • ||i is the Schatten 1-norm defined in 
Appendix [X] 

It should be noted that an alternate measure of distance between density ma- 
trices to that provided by the Schatten 1-norm is given by the fidelity function |23j . 
One could write an alternate form of the QCC in terms of the fidelity that would be 
essentially equivalent to the form of the QCC given in 1(12(1 above. Using an obvi- 
ous notation to denote the QCC written with each of these definitions of distance, 
the two forms are related as follows. Given a quartet {P, U, M. enc , -Mdec}, one can 
show that if QCC Schatten (P, U, A^eno -A4dec> a) is satisfied, then the fidelity-based 
version of QCC given by QCC fidelit (P, U, M e nc, -Mdec, a') is also satisfied, where 
a' = a' (a) is a function of a. The form of QCC based on the Schatten 1-norm given 
above in 1)12(1 is more convenient for our purposes for a number of mathematical 
reasons, including the fact that that the fidelity, as such, is not a proper norm. For 
instance, the statement and proof of the Encoding No-Go Theorem carried out in 
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fJSlbclow are more conveniently presented making use of the form of the QCC based 
on the Schatten norm. 

Note that Alone and Aldoc do not represent physical operations: all physical 
operations are carried out by the quantum component P. If the proper distinction 
between these maps and P is not observed, one could include the entire computa- 
tion in the definition of the maps, with the absurd conclusion that any quantum 
computation could be performed in the absence of any real hardware. 

2.3.2. Some implications of the QCC. The QCC is a remarkably powerful ex- 
pression. It constitutes a kind of "master expression" for physical quantum com- 
putation. The inequality l|12|l concisely incorporates a complete specification of 
the full dissipative, decohering dynamics of the actual, practical device used as the 
quantum computing machine, a specification of the ideally-defined quantum com- 
putation intended to be performed by the machine, and a quantitative criterion for 
the accuracy with which the computation must be executed given the inevitability 
of residual errors surviving even after error correction has been applied. 

Making use of the QCC, one can state and prove (we do this is in JO the next 
section of the paper) a fundamental and powerful theorem in the subject of quan- 
tum computing, the Encoding No-Go Theorem. This no-go theorem gives a pre- 
cise mathematical statement of the conditions under which fault-tolerant quantum 
computation becomes impossible in the presence of dissipation and/or decoherence. 
Apart from its formal significance, the theorem can be used to compare different 
proposed physical approaches to actually building a quantum computing machine, 
with the no-go condition furnishing a criterion for the practical engineering viability 
of various choices. 

As a further indication of the power and general applicability of the QCC, 
we show that one can apply it to the known, seemingly distinct "paradigms" for 
quantum computing, based on (1) the use of quantum circuits built up out of 
quantum gates (the circuit-based paradigm, or "two-way" quantum computing), 
(2) the use of graph states or cluster states (the graph state-based paradigm, or 
"one-way" quantum computing) and (3) the use of specially chosen Hamiltonians 
describing adiabatic dynamics (the adiabatic quantum computer paradigm). The 
QCC allows one to show that these apparently different definitions of a quantum 
computer are in fact manifestations of the same underlying formulation: there is 
only one paradigm for quantum computers. The application of the QCC to different 
quantum computing paradigms is discussed in 3SI 

The encoding-decoding pair Af C nc and A^dec that appear in the QCC are de- 
fined quite generally as completely positive trace-preserving maps. This formula- 
tion is sufficiently general to encompass all possible encodings associated with stan- 
dard quantum error correction (QECC) techniques, dccohcrence-free subspaces and 
noiseless subsystems. More generally still, we show below in t!2.3.3l that the recently 
discovered approach known as "operator quantum error correction" (OQEC) is in 
fact a special case of our more general QCC formulation. In addition the QCC 
can be used to extend OQEC to what we will refer to as "operator quantum fault- 
tolerance" (OQFT). In particular, we show in ^ below how the QCC allows one 
to derive error thresholds for fault tolerance in a completely general context. 
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Another significant consequence of the QCC is that it resolves the apparent 
paradox that the quantum computations we wish to perform are defined by re- 
versible operators, but the actual devices that we must use to execute the compu- 
tations are necessarily described by irreversible maps. We note that this is remi- 
niscent of the fluctuation-dissipation theorem, which relates irreversible dynamics 
to equilibrium properties in a large class of physical systems. Inspection of l|12(l 
reveals that the paradox is resolved through the transformations provided by the 
encoding and decoding maps associated to the QCC. Roughly speaking, reversible 
behavior of the actual device is enabled only on the code subspace defined by A4 e nc 
and Mdcc- 

Note that if one describes quantum computing solely in terms of the unitary 
transformations that define the computations, it is not too surprising that the 
resulting computational model is in some way "powerful." After all, unitary trans- 
formations on finite dimensional spaces include such powerful operations as the 
discrete Fourier transform which are known to play an important role in num- 
ber theoretic problems. Rather than regarding the power of the ideally-defined 
quantum computation as the remarkable thing, the truly remarkable thing would 
be the construction of an inherently irreversible device that actually implements 
the reversible, unitary map to a specified level of accuracy. It is this possibility 
that our QCC expresses and makes mathematically precise. The QCC can thus be 
used to formulate and prove assertions about the physical solvability of particular 
computational problems. 

2.3.3. Operator quantum error correction (OQEC) as a special case of QCC. 
The recently developed theory of operator quantum error correction [18) , |19j uni- 
fies many apparently disparate approaches to the practical problem of dealing with 
errors in quantum information. Among these approaches are quantum error cor- 
rection, decoherence free subspaces and noiseless subsystems. Here we show that 
OQEC is in fact a special case of the general statement of the QCC. In this section 
we demonstrate that the entire formalism of OQEC can be obtained from QCC by 
choosing M. enc and A'idec as described below, and by setting setting U = I and 
a = in the QCC, so that the reduction QCC — * OQEC is given by: 

(13) QCC(P,U,M cnc ,M dcc ,a) -> QCC(P,I,M cn c,M doc ,0) . 

In the scheme of OQEC, the techniques of (standard) quantum error correction, 
decoherence free subspaces and noiseless subsystems are subsumed under the unified 
concept of "correctability." This is defined formally as follows. Let H comp be a 
Hilbert space with a decomposition H comp = (H A <g> H B ^j ©A', where H A , H B and 
K are discussed in the following paragraph. We identify this as the computational 
space -ffcomp of this paper because, as we shall see below, it describes the space 
on which the real physical processes of error and recovery operate. Let 6 be the 
semigroup given by 
(14) 

6 = {<tE L(H comp ) : 3cr A G L(H A ), 3<r B G h(H B ), such that a = a A ® a B } , 

where L(-) is the space of bounded operators 6 on the appropriate Hilbert space, with 
the operator norm, and let £ and 1Z be completely-positive, trace preserving maps 



6 In i|2.3.c{l of this paper wc assume (following that all Hilbert spaces are finite- 

dimensional. Although not discussed in 1191 . it is important to note that if one makes this 
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on L(i7 com p) corresponding to error processes and recovery processes, respectively. 
We say that 6 is correctable for £ iff 

(f 5) Vcr G 6, (Tr^ oV 6 oTZoS)a = Ti A a , 



where Ve effectively projects density matrices onto the H £g> H B subspace of 

Physically speaking, H B is the Hilbert space which carries the information to 
be protected from errors, the "noiseless subsystem," while H A is the Hilbert space 
on which the errors are permitted to operate freely, the "noisy subsystem." K is 
the orthogonal complement of H A ® H B in H corav and is simply projected out 
by Ve- Thus, in 1)1 Tr^cr represents the quantum information which is to be 
protected. The left side of (|15)l describes the effect of allowing errors to operate 
on the full state a, and then applying recovery procedures. (The projection and 
the trace simply extract the state of the noiseless subsystem.) According to (fT~5f> . 
correctability thus means that the recovery procedure does in fact recover the state 
of the noiseless subsystem, Tr act, without error. 

As we now show, the definition of correctability is actually a special case of the 
QCC. Note first that correctability, as defined above, applies to a quantum channel 
as opposed to a quantum computer. It describes the transportation of a quantum 
state in a noisy environment as opposed to the "processing" of a quantum state so 
as to implement a quantum computation. In order to make the connection with 
the QCC, we may therefore think of the quantum channel as a quantum computer 
that implements the identity operation: 

(16) U = / Hlogical 

In the definition of correctability, the part of the state containing the informa- 
tion of interest is recovered without error. This corresponds to taking a = in the 
QCC. Note that, as discussed above, this is not practically achievable for quantum 
computers. This is less obviously an issue for the theory of operator error correction 
as currently formulated |18j . |19| . since that theory addresses a relatively circum- 
scribed set of circumstances in which one is concerned with transporting quantum 
states rather than implementing a quantum computation. In particular, the er- 
ror process £, which is specified in advance, only operates once on the quantum 
state being transported. Error processes associated with the constituent parts of 
quantum computers operate each time the constituent part operates on the the 
quantum state being processed. The problem of fault tolerant quantum computa- 
tion is inherently more complex than the problem of error correction/protection for 
a quantum channel. We will discuss this more fully in ^Hl below. 

Setting U — I and a = in the quantum computer condition, we obtain 

(17) HMd BB °P°M mD )p-p\\ 1 = 0, 
or 

(18) (Mdec ° P ° Menc) P = P , 

where p £ L(Hi ogical ). 7 

assumption, there is then no need to distinguish between bounded operators in L(-) and trace- 
class operators in T(-). This distinction is important in the other sections of our paper where we 
allow Hilbert spaces of infinite as well as finite dimensionality. 

^As noted in the previous footnote, we are assuming in H2.3.3I that -ffi og ; ca i is finite- 
dimensional. 
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We now proceed to define the encoding map Alone- For this purpose we define 
a map W cnc ■ L(-ffiogicai) — * L(i7 B ) that encodes the logical quantum state p in a 
state a B of the noiseless subsytem. We further define a map W ad j(a A ) : L(H B ) — > 
li(H A <g> H B ) that adjoins an arbitrary state a A of the noisy subsystem to the state 
a B , that is 

(19) W adj (a A ) : a B ^ a = a A <g> a B 

We then define the full encoding map that appears in the QCC as follows: 

(20) M enc = W adj (<J A )oW cnc 

The map P characterizes the dynamics of the physical computer, which in this 
case is just the (noisy) channel followed by the recovery procedure: 

(21) P = UoS 

We note that the current formulation of the theory of operator quantum error cor- 
rection implicitly assumes that the recovery process 1Z can be implemented without 
error, even though this requires, in general, that coherent operations be performed 
on entangled quantum states. Once again, the emphasis on error correction alone 
avoids the more difficult issue of achieving true fault tolerance. 
Finally we define the decoding map as: 

(22) M dcc = W e ~l o Tr A ° V 6 , 

which extracts the state of the noiseless subsystem and decodes it to obtain a state 
in the logical space L(ffi gi C ai)- 

With the above definitions, the QCC becomes 

(23) (W~l o Tr A o V e ° TZ o 8 ° (a A ) o W cnc ) p = p. 
Applying W cnc to both sides, and recalling that W cnc p = <j b = Tr^er, we obtain 

(24) (Tt a oV e oTZo£o W adj (a A ) o W cnc ) p = Tr A a . 

Noting that Wa d j(<j A ) o W cnc p = cr, we obtain the correctability condition of (18] . 

m 

(25) (Tr A oV & oTZoE)a = Ti A a . 

In summary, we have shown that the formalism of operator quantum error 
correction actually arises as a special case of an underlying definition of physical 
quantum computation given by the QCC. In addition, examination of operator 
quantum error correction (OQEC) from the perspective of the QCC reveals limi- 
tations and restrictions implicit to the formalism of operator quantum error cor- 
rection, and inherited from the previously known, standard approaches to error 
correction (QECC). These limitations render direct application of either OQEC or 
QECC to questions of fault tolerance somewhat problematic, whereas the QCC 
approach is more immediately applicable. Thus, the QCC enables one to generalize 
OQEC to operator fault tolerance (OQFT). 
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3. The Encoding No-Go Theorem 



3.1. Introduction. Armed with the QCC, in this section we state and prove 
a theorem of crucial importance in the theory of physical quantum computation. 
This is the Encoding No-Go Theorem, which gives a precise mathematical statement 
of the conditions under which fault-tolerant quantum computation becomes impos- 
sible in the presence of damping. Damping, for which we provide a mathematically 
rigorous definition below, includes the effects of dissipation and decoherence. The 
No-Go theorem for a completely positive trace-preserving map P corresponding to a 
putative quantum computing device then asserts that, in the presence of sufficient 
damping, the quantum computer condition QCC(P, U, A4 enc , .Mdcc, a) (cf Q12|ll 
cannot be satisfied for any encoding-decoding pair, unless -ffi og icai has dimension 1: 
there is then effectively no quantum computer. (In the case that dim illogical = 1 we 
are of course unable to define a meaningful quantum computation at all.) As part 
of the No-Go Theorem we explicitly calculate a universal critical damping value for 
fault-tolerant quantum computation. We precisely state and prove this theorem in 
the remainder of £0 

3.2. Encoding and Decoding Maps. An encoding-decoding pair are com- 
pletely positive, trace-preserving maps 



Encoding-decoding pairs provide the link between a completely positive map P : 
T(-ffcomp) —* T(i7 comp ) corresponding to a physical device and a unitary operator 
U : i?iogicai — * -/^logical corresponding to a quantum computation. We will place no 
further restrictions on encoding-decoding pairs. Now, suppose .Menc-Mdcc is an 
encoding-decoding pair. Then the adjoint maps .M dcc , M\ nc are unit preserving 
completely positive maps. (Adjoint maps of completely positive maps are defined 
in (|T32|> in Appendix |B]) 

3.3. Damping and 7-damping. Dissipative and decohering effects are con- 
sequences of damping. To begin the development of the No-Go Theorem, we intro- 
duce a mathematically precise definition of the damping of a quantum mechanical 
system, to which we refer as 7-damping. 

Definition 3.1. Let H be a Hilbert space. A completely positive trace- 
preserving map P : T(H) — > T(iJ) is ^-damped iff there is an abelian von Neumann 
algebra 21 C h(H) such that for all T 6 L (//"), there is an St G 21 such that 



where the operator norm || • jjoo is defined in Appendix lAl Note that larger values 
of 7 correspond to less damping of the system. 

To make contact with the physically intuitive notion of damping, we apply 
this definition to the example of the simple harmonic oscillator subject to phase 
damping. For such an harmonic oscillator, the ijth element of the density matrix, 
Pij = (i\p\j), decays exponentially as e~ K ^ l ~^ , where we are working in the basis of 
energy eigenstates identified by the labels i and j. The quantity k is characteristic of 
the specific oscillator and its coupling to the environment. The completely positive 
map P transforms the initial, general density matrix for the system into a density 
matrix with exponentially-decaying off-diagonal elements. Under the influence of 



(26) 




(27) 



\\P\T)-S T U<l\\T\ 
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damping, as the off-diagonal states of the oscillator decay and approach zero, the 
density matrix converges to a diagonal density matrix in the {ij} basis specified 
above. This is true for all initial configurations of the oscillator, and thus the 
damping process tends to an abclian set of final configurations. (The damping 
parameter k that characterizes the decay of the off-diagonal elements of the density 
matrix is related to the quantity 7 that appears in l|27|l. As noted above, larger 
values of 7 correspond to less clamping and hence smaller values of n.) 

3.4. No-Go Theorem for Encodings. We now state the main result of fJ3 
the Encoding No-Go Theorem. 

Theorem 3.2 (The Encoding No-Go Theorem). 

Suppose that QCC(P,U,M enc ,M de c,a) holds. If P : T(H comp ) -> T(i7 comp ) is 
^/-damped and 27 + a < a/2/4, then /fi og icai has dimension 1. 

To prove the Encoding No-Go Theorem, we first derive in H3.5l a number of general 
mathematical results on completely positive maps. We then apply these specifically 
to obtain a proof of the Encoding No-Go Theorem in H3. 61 below. 

3.5. Completely Positive Maps with Abelian Factorizations. We begin 
with the following lemma, which furnishes a superoperator version of the definition 
of 7-damping. For any abelian von Neumann algebra 21 C L(ff) there is a unit- 
preserving completely positive projection operator E<a : — » 21. This fact is 
elementary, but also follows from injectivity [7] of such algebras, in the case H is 
separable. 

Lemma 3.3. If P is "/-damped, 21 is as given in Dehnition Vd. 1\ and Egi : L(i?) — > 
21 is a unit-preserving completely positive projection mapping E^ : L(_ff) — > 21, then 

(28) WP'iT) - EaF^TJIloo < 27||r|| 00 . 
Proof. Note that Ea is a linear mapping of norm < 1 and thus 

\\P\T) - E<nP t (T)\\ oc < \\P\T) - St\\oo + \\S t - E^COIU 

(29) < -yllriloo + ||E a( ST-E 2[ P t (T)|| 00 

<27l|Tl|oo, 

as claimed. □ 

In the remainder of this section we prove an important technical result used in 
the proof of the no-go theorem, namely that completely positive maps F : L(i?) — > 
Jj(H), which factor through completely positive maps into abelian von Neumann 
algebras, cannot be used to approximate unitary operators U on H if dim(iJ) > 2. 
More precisely, we will show that for any (3 < \/2/A, there is at least one non-zero 
operator T for which 

(30) ||L/ t TC/-P(T)|| 00 >/3 llTHoo. 

We first define what it means for a completely positive map to factor through an 
abclian von Neumann algebra: 

Definition 3.4. Let A be a C*-algebra with a multiplicative unit. A unit 
preserving completely positive map F : A — > A has an abelian factorization (or 
briefly is abclian factorizablc) iff there is an abclian von Neumann algebra 2} and 
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unital completely positive maps Q : A — ► 03, R : 03 — > A such that F is the 
composition Ro Q. 

If F : L(i7) — > li(H) is abelian factorizable, then it follows from the definition 
that for any unit-preserving completely positive maps Q, R, the completely positive 
map Q o F o R is abelian factorizable. We begin by providing a characterization of 
abelian factorizable maps. 

Proposition 3.5. Let A be an arbitrary C* -algebra with multiplicative unit. 
A unit preserving completely postive map F : A — > A factors through a finite- 
dimensional abelian von-Neumann algebra iff there are unit preserving positive lin- 
ear functionals p\ , . . . , p m G A* and positive elements G\, . . . , G m E A of norm 
< 1, such that y^ ? -_ 1 Gi = I and 

m 

(31) F(T) = J2pi( T ) G i VTeA 

i=l 

PROOF. We first note that positive maps from an abelian C*-algcbra or into 
an abelian C*-algebra are automatically completely positive (see |29|.[(j] ). Thus 
the map F given by Equation (|31|) is completely positive. Let 05, Q : A — > 25, 
R : 05 — > A as in Definition 13.41 but with 05 finite dimensional and let E\ , . . . , E m 
be the minimal non-zero projections of 03. Then 



!=1 



(32) F{T) = R(J2E t Q(T)E t ) = R^Pii^E, ) = ^ Pl {T)R(E l )- 
Letting G t = R{E l ). 

m 

(33) y. g 1 = H r ^) = r [H e A = L 



This completes the proof. □ 

In general, unit-preserving completely positive maps with arbitrary abelian 
factorizations can be approximated by maps of the form (|31|l . 

Proposition 3.6. If a completely positive map F : A — > A has an abelian 
factorization, then there is a generalized sequence of maps {F k } k ^k each having 
the form F K (T) = 2j=i Pi{T)Gi which converges to F in the point-norm topology, 
that is for each T G A, F K (T) — > F(T) in the norm of A. 

Proof. Let 05, Q : A -> 03, R : 03 -> ^ as in Definition EH Given 
Ti, . . . , T m G A, let 05o be a finite dimensional abelian von-Neumann subalgebra of 
05 and E a linear projection operator 8 03 — > 05o such that 

(34) ||Q(T i )-E(Q(T i ))|| 00 < e. 
Since Q is contractive, it follows that 

(35) \\F(Ti) - i? o E o Q(Ti) Hoc < e. 

Now apply the preceding result. □ 



8 These operators are sometimes referred to as conditional expectations. 
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Assume H is a finite dimensional Hilbert space. We will consider trace func- 
tionals on two distinct spaces of operators: one on the space Jj(H), which we denote 
by tr, and the other on the space L(L(i?)) of linear mappings L(i?) — > L(if), which 
we denote tr oper . We will prove that abelian factorizable completely positive maps 
cannot approximate the identity map on L,(H). To do this we will show that the 
trace functional tr opcr separates, in a sense to be made precise in the next para- 
graph, the identity operator on L(iJ) from abelian factorizable completely positive 
F. Note that tr opcr (7 L(H) ) = dim 2 (ff). 

Lemma 3.7. Suppose the Hilbert space H has finite dimension n. For any 
unit-preserving abelian- factorizable completely positive map F : L(_ff) — > L,(H), 

(36) tr oper (F) < n. 

Proof. It suffices to prove this for F which have the form (|31() . Referring to 
that representation, each positive functional pi can be represented by a non-negative 
operator Si as follows: 

(37) Pi (T)=tv(TSi), 

Since Pi(I) = 1, Si also has unit trace and in particular, Si < I. Moreover 
EtiG l =Ih- Now 

m 

tr oper (P) =£tr(S i G 1 -) 

m 

= ^(Gl /2 S z Gl /2 ) 

(38) Z 1 

<$>(GW) 

i=l 
m 

= 5^tr(G0=tr(I H )=n. 

»=i 

□ 

The previous lemma gives us a lower bound on the trace of I — F for F abelian 
factorizable: 

(39) tr opcr (7 - F) >n 2 -n. 

We can use the above lower bound on the trace of /— F to obtain a lower bound 
on the norm of the operator I — F : L(iJ) — > ~L(H), where we consider L(ff) with 
the Schatten 2-norm || • 1 1 2 , defined in Appendix 1X1 The Schattcn 2-norm is the 
norm that arises from the trace inner product on L(77), also known as the Hilbcrt- 
Schmidt norm. We denote the corresponding operator norm on L(iJ) — > L(iT) by 

II -H2-2. 

If F is self-adjoint as an operator on the space L(i?) with the trace inner 
product, then from Lemma \'S. 71 we immediately obtain the bound 

(40) ||7-F|| 2 -2>l-i. 

n 

In fact, the lower bound l|4U|l is true for general unit-preserving abelian factorizable 
completely positive maps F. To see this, write F = F^ t + iFy m where both Ftn e , 
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Fy m are self- adjoint operators (not necessarily completely positive, however). Now 

(41) tr oper (/ - F %t ) = tr opcr (I - F) > n 2 - n. 
Therefore 

(42) ||7 - F\\ 2 ^ 2 >\\I- F mt \\ 2 ^ 2 > 1 - 1/n. 

In the discussion that follows, we need to consider another norm on the space 

L(iJ) — > L(_ff) in addition to the norm || ■ 1 1 2 >2 just considered. The new norm, 

which we denote || ■ ||oo^oo; is also an operator norm on L(i7) — > L(iJ), but relative 

to the || ■ |oo norm on L(H). The || • ||oo^oo norm is different from the || • 1 1 2 ^2 norm, 

but for finite dimensional spaces H the two norms are equivalent, which means that 
each norm is bounded relative to the other. To obtain the bounding constants, note 
that if T e L(H), 

(43) imu < TO < v^imioc. 

From this it follows that 

(44) 4= H F ll2-2 < II^Hoo^co < V^!l^l|2^2. 

which is the desired relative bound. In 1)40(1 . if we substitute F by I — F, we 
immediately obtain the following proposition: 

Proposition 3.8. Suppose H is a Hilbert space of finite dimension n. If a 
unit-preserving completely positive map F : L(77) — > L(i7) has the form (|3 1 1> - then 

(45) \\I-F\\ 00 ^ 00 >n- 1 ' a (l-l/n). 

To prove the crucial result for the No-Go Theorem, we only use case n = 2 
of (35) . 

Theorem 3.9. Suppose H is of dimension > 2. If F is a unit-preserving 
completely positive map on L(if) with an abelian factorization and [3 < V2/4, then 

(46) ||tftTl7-.F(T)|| 00 >/3 \\T\\ X 

for at least one T € L*(H). 

If H is finite dimensional, we can take (3 = y/2/4. 

Proof. Replacing F by the completely positive map T 1— > UF(T)U\ we can 
assume without loss of generality that U = I. To prove this, we will show that the 
assertion that 

(47) HT-^TJHoo <(3 \\T\U VTeL(H), 
leads to a contradiction. However, (|47|l implies 

(48) ||J - i^Hoo^oo < /3. 

We now reduce the proof to the case H has dimension 2, by considering a Hilbert 
space K of dimension 2 and completely positive unit-preserving mappings Q : 
L(K) -> L(H) and R : L(H) -> L{K) such that RoQ is the identity map on L(K). 
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If Ih can be approximated to within f} of the abelian factorizable F, then RFQ is 
also abelian factorizable and 



(49) \\I K - RFQ\\ = \\RI H Q - RFQ\\ < \\ I H - F\\ < [3 < V2/4. 

However, in this contradicts Proposition ^. 81 

In the finite dimensional case, the norm is actually achieved so that in l|48|) the 
< sign can be replaced by < and so we can take (5 < as claimed. 9 □ 

3.6. Proof of Encoding No-Go Theorem. 

3.6.1. Two Lemmas. 

Lemma 3.10. Suppose that QCC(P, U,M cnc , Mdca ") holds. IfP: T(iJ comp ) - 
T(-ffcomp) is ^-damped and 21 and E<a are as identified in Lemma Vj. 31 then for every 

gical ) } 

(50) WMl^EvP'M^T) - rfTUWoa < (2 7 + o)||T||oo. 

PROOF. The QCC(P, U, M cnc , Mdcc, a) implies that for every p £ T(iJi ogical ) 
with ||p||i < 1 and self-adjoint T £ L(iJi ogica i), 
(51) 



<M Acc {P(M cac (p)))-UpU^\T = tr 

< <A\ T 

It follows that for every T £ L(iJi ogica i), 

(52) llKnc(^WdecCO)) - f/^C/Hoo < a||T|| 

21 is commutative and by Lemma 13.31 for each T £ L(7Ji og i ca i) 

(53) HP^T) - EaP^T)!^ < 2-ylirHoo. 

Thus using the fact that Mdec & n d M enc have norm < 1, for every T £ L(i7i og i ca i), 

WMl^P'M^Tj-U^TUWoe 

< WMl^P'M^iT) - M^P'M^T)^ 
^WM^P'M^m-^TUU 

< (2 7 + a)||T|| 00 . 

□ 



p\Ml nc {P\M\ ec {T)))-U^TU 



| oo • 



(54) 



Lemma 3.11. If -ffi og icai is of dimension > 2, and P, U , A4 cnc , AA& CC and 
Ea are the same as in Lemma \3.KA and if (3 < v^/4, then for some non-zero 
T £ L(i?i og i ca i) ; 

(55) ||ML c Ea^decCn - rtTUU > (3\\T\\oo 

PROOF. The unit-preserving completely positive map R = M t cnc E%P t M dcc 
factors through the abelian von Neumann algebra 21; this follows from the presence 
of the projection operator Ea in the expression for R. Then l|55|l follows from 
Theorem ESI □ 



^We note that it is possible to derive explicit results as well for the case n = 3. In this 
case one can show that one obtains a larger numerical bound than V2/4, which applies when 
dim(ff logical ) > 3. 
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3.6.2. Statement of Proof of Encoding No- Go Theorem. 

PROOF. By the hypotheses of the Encoding No-Go Theorem (Theorem 13.20 . 
the quantity 27 + a < \/2/4. Choose f3 such that 2j + a < (i < v2/4. From Lemma 
I3.1UI it follows that for every T £ L(//i og icai), 

(56) WMt^E^Ml^T) - UtTUWn < (2 7 + a)||T|U. 

On the other hand by (|55|l in Lemma \'S. Ill there is a non-zero T such that 

(57) IIMLcEaP^ecCn - U^TUU > /3HTIU. 

Since ||T||oo > 0, equations l|56[) and (|57|l imply (3 < 27 + a, which contradicts the 
choice of (3. □ 

3.7. Interpretation of Encoding No-Go Theorem. The Encoding No-Go 
Theorem is an extremely powerful application of the QCC. With this theorem, one 
can calculate the amount of damping for which fault-tolerant quantum computation 
becomes impossible. When the amount of damping exceeds the critical amount, 
which means that the value of 2 7 becomes less than the critical value 27 cr iti ca i = 
-\/2/4 — a, we find that the only solutions to the QCC are those for which dim 
-ff logical = 1. As noted above, in this case no meaningful quantum computation is 
possible, since not even one quantum bit can be accomodated. 

In order to analyze the constraints on solutions to QCC(P, U, M C nc 7 A4dcc, a.) 
implied by the No-Go Theorem, we regard the pair {U, a} as given, since both 
the desired quantum computation U , and the maximum acceptable implementa- 
tion inaccuracy a, are prescribed. We assume that U is defined on a Hilbert space 
-^logical such that dim (i7i gi ca i) > 2, to allow meaningful quantum computation. 
For practical applications, one then seeks to determine triples {P, A4 enc , Mdec} 
that satisfy QCC(P, U, A4 enc , -Mdea a). The No-Go Theorem, on the other hand, 
identifies conditions in which QCC(P,U,M enc ,Mdcc,o:) is not satisfied. Thus 
the No-Go Theorem provides a useful calculational tool for eliminating prospec- 
tive quantum computer implementations that are guaranteed to fail. With the 
criterion provided by the No-Go Theorem, we can bound the space of solutions 
to QCC(P, U, Monc, Mdcc, ct). This allows exploration of trade-offs amongst the 
members of the triple {P, M e nc, A^dec}, with which one may construct general- 
ized "phase diagrams" that indicate boundaries between potentially acceptable and 
definitely unacceptable values for P, A'fenc and A^dec- 

4. Quantum Components 

4.1. Introduction. As defined in Sectional we refer to a physically realizable 
device intended to implement a quantum computation as a quantum component. 
Mathematically, a quantum component is represented by a completely positive, 
trace preserving map, P. The detailed form of P, as an explicit function, is dictated 
by underlying equations of motion. 

For a closed physical system, the quantum mechanical dynamics of the system 
are given by the Schrodingcr equation associated to a particular Hamiltonian Tl. 
However, for the general problem of a practical quantum computer, we must analyze 
realistic quantum components that interact with their environment, dissipate heat 
and exhibit decoherence. We must thus utilize a formulation that yields equations 
of motion appropriate to open quantum mechanical systems. 
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The connection to the Quantum Computer Condition presented in Section El is 
made by starting with the appropriate equations of motion that describe the dynam- 
ical evolution of a realistic quantum component interacting with its environment. 
One proceeds by solving the appropriate equations of motion. The explicit solution 
thus obtained provides the time-dependence of the quantum mechanical state of 
the open system. In principle, this allows us to deduce the explicit functional form 
of P. 

Formulations of equations of motion for quantum components that interact with 
complex environments comprised of many degrees-of-freedom necessarily involve 
approximations of one sort or another, and there is not in general a unique choice. 
In this section we illustrate the general approach by making use of a Lindblad- 
type equation (made more precise below) to describe how one obtains the quantum 
component P that appears in the Quantum Computer Condition. 



4.2. Dynamical Equations of Motion. 

4.2.1. Time- dependent generalization of the Lindblad equation. The state of 
a quantum mechanical system defined on a Hilbert space H can be modeled by 
a density operator p(t) whose time dependence obeys a linear (time dependent) 
equation 

(58) lL p{t) = A(t)p(t), 

where A(t) is an operator acting on T(iJ), the Banach space of trace class operators 
on H. For a closed system we have the Schrodinger equation, and the right-hand- 
side of H58|) is given by 

(59) A(t)p = ~[H(t),p], 

where the Hamiltonian Ti{t) is a self-adjoint operator which may depend on the 
parameter t. 

However, it would be impractical to describe realistic quantum components 
using (|59f) . since writing down the detailed Hamiltonian operator to account for all 
of the degrees-of-frccdom comprising the quantum component and its environment 
would be intractable. 

Motivated by the work of Lindblad |21j , we make use instead of the following 
expression for the action of A(t) on p: 

(60) A(t) P = - l -[n(t),p] +J2 k-(*)p4(*) - \{ L ]V) L Mp) 

3 

where as above the Tl(t) is a self-adjoint operator which may depend on the pa- 
rameter t, and the Lj are operators that describe effects arising from interaction 
with the environment, such as dissipation and decoherence. These operators are 
generalizations of the Lindblad operators of |21j ; unlike the treatment given by 
Lindblad in [21] . however, which considered only time- independent Hamiltonians 
7i and time- independent dissipative perturbations Lj, with all operators bounded, 
we will allow unbounded and time-dependent Ti.(t) and time-dependent (but still 
bounded) dissipative perturbations Lj(t). Our treatment generalizes that given 
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21] . and we will refer to the equation of evolution 

L (t)pL]{t)-\{L]{t)L 3 {t),p} 



(61) j t p& = -\m),p\+Y. 

j 

as the time- dependent generalization of the Lindblad equation. As an example of 
how one may approximate the dynamics of a quantum component interacting with 
a complex environment, this equation provides a starting point to the derivation of 
P used in the Quantum Computer Condition. For this purpose we need to consider 
solutions to equations of the type given in (|61|l known as "fundamental solutions." 

4.3. Fundamental Solutions. The concept of a fundamental solution asso- 
ciated to an evolution equation (see |30j . §4.4) of the general form l|58|) . such as 
(16111 in particular, where p is a function with values in the Banach space T(H) and 
A(t) is a one-parameter family of (possibly unbounded) linear operators on T(H), 
will play an important role in this paper. A fundamental solution associated to an 
equation of motion is a solution of an operator version of the original equation. We 
show below how, given a fundamental solution, one obtains the completely positive, 
trace preserving map P that appears in the Quantum Computer Condition. 

Fundamental solutions are given by a family {Pt,s}t>s>o of bounded operators 
on T(H) indexed on pairs of real numbers t and s satisfying equations 

(62) j f Pt,s = A(t)P t>s 
and 

(63) P s . s = I. 

The intended interpretation of Pt iS is that if the system is in state p at time s, then 
the system will be in state Pt }S • P at later time t, that is 

(64) p(t) = P t . s ■ p(s). 

4.3.1. Existence and Positivity Properties of Fundamental Solutions. The prob- 
lem of the existence and uniqueness of fundamental solutions is a central one in the 
mathematical theory of evolution equations. In addition to addressing the ques- 
tion of the existence of fundamental solutions, it is important for our analysis to 
determine the positivity properties of fundamental solutions. This is because, as 
we shall see, the complete positivity of the map P that explicitly appears in the 
Quantum Computer Condition is inherited from the complete positivity of the set 
of fundamental solutions {Pt iS }. Positivity is important in ensuring that the set 
{Pt iS }, as well as P, carry density matrices to density matrices. 

For equations of motion associated to finite dimensional systems, existence and 
uniqueness of solutions follows from the Lipschitz theorem on ordinary differen- 
tial equations. Lindblad's analysis extended to infinite-dimensional (but bounded) 
systems, so, more generally, Lindblad characterized the infinitesimal generator of 
a norm continuous completely positive semigroup |21j . which corresponds to the 
case of <|58[1 in which A(t) is constant and norm bounded (but possibly infinite- 
dimensional), i.e., for the standard, time-independent Lindblad equation. 

In order to generalize Lindblad's analysis to allow us to study infinite-dimensional, 
unbounded, time-dependent quantum systems, we will need to consider general evo- 
lution equations in which A(t) may be unbounded as well as time-dependent, 
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for the infinite-dimensional case. The general theory of such equations was devel- 
oped by Kato, Yosida and others in the 1950's. We will rely on results of Kato [14] 
which pertain to both existence of solutions and positivity properties, and on The- 
orem 4.4.1 of |30| which pertains to existence of solutions; moreover, we will use a 
constructive form of the theorem (which follows from an examination of the proof) 
which expresses the fundamental solution as a limit of a product of exponentials. 
The basic assumption of the approach is that if the A(t) are generators with suffi- 
ciently smooth variation, then fundamental solutions exist. 

Our system comprised of a quantum component interacting with its environ- 
ment, described by l|61(l . consists of a time-dependent Hamiltonian Ti(t) character- 
ized by a time-dependent perturbation V(t) of a (possibly unbounded) self-adjoint 
operator Ho so that we have H(t) = Ho + V(t). In order to apply Kato's theory for 
time dependent evolution equations, we will assume among other things that the 
perturbation V(t) does not change the domain of Ho- For the necessary background 
see f |14| .[30 ] l. We now state a proposition that asserts the existence and complete 
positivity of fundamental solutions to operator versions of the time-dependent gen- 
eralization of the Lindblad equation given in (|61|) : 

Proposition 4.1. For suitably regular time-varying potentials V{t) and dissi- 
pation operators Lj(t), there exists a strongly continuous completely positive oper- 
ator P t , s which is a fundamental solution to the time- dependent generalization of 
the Lindblad equation. 

The exact statement of the conditions for the above in the form of a theorem, 
and proof, are given in Appendix IBI in ffl3.3l and ffl3.4l 

4.4. Construction of Quantum Components. Having established the ex- 
istence and complete positivity of fundamental solutions to the operator form of the 
underlying equation of motion for our system (comprised of the quantum compo- 
nent interacting with its environment), it is straightforward to obtain an expression 
for the completely positive, trace preserving map P, characterizing the quantum 
component, that explicitly appears in the QCC. This is simply obtained by noting 
(cf jnUl) that the time-evolution of the state p(t) between an initial fixed time s 
(corresponding to the start of the quantum computation) and a final fixed time t 
(corresponding to the end of the quantum computation) is fully specified by the 
fundamental solution defined with respect to those time values, P^g, so that we 
have 

(65) p(t) = P Vs ■ p(s), 

and hence the completely positive, trace-preserving map P that appears in the 
QCC is given by the equivalence 

(66) P = Pu ■ 

5. Unified Treatment of Quantum Computing Paradigms 

5.1. Introduction. In this section we show that the QCC provides a uni- 
fying framework in which to describe on the same footing the currently-known 
"paradigms" for quantum computation, including the circuit-based paradigm, the 
graph state-based paradigm, the adiabatic quantum computer paradigm. The QCC 
subsumes all of these into a single, unifying paradigm for quantum computing. 
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5.2. Circuit-based paradigm. In this section we describe the specification 
of quantum components based on the "circuit-based" paradigm of quantum com- 
putation. We proceed as follows: 

(1) For purposes of clarity, we begin in l5.2.1l bv working in an idealization in which 
there is no noise present. For this idealized case we obtain the general form of 
the completely positive map P characterizing the quantum component. We then 
apply the result (for the noiseless idealization) to several specific realizations of 
the circuit-based paradigm. These include qubit-based quantum computers (these 
utilize states, the operators for which have a discrete eigenspectrum, i.e., qubits in 
the case of 2-level systems), quantum continuous variable-based quantum comput- 
ers (these utilize states, the operators for which have a continuous eigenspectrum, 
such as coherent states), and liquid state NMR-based quantum computers. 

(2) Having obtained the general form for P in the noiseless case, for each of the 
three above-mentioned realizations of the circuit-based paradigm, we then explain 
in 15.2.21 how to modify the analysis, in a way appropriate to all choices of circuit 
realization, so as to account for the effects of noise. 

5.2.1. Idealized Circuits in the Absence of Decoherence and Dissipation. A quan- 
tum circuit is described by a set G of gates operating in some specified order on 
elements of a set R of objects. The quantum states of these objects constitute the 
information which is "processed" by the gates of the quantum circuit. We associate 
with each object i of R a Hilbcrt space that describes the possible states of 
that particular object. The Hilbert space for the full set of objects on which the 
circuit operates is then 

(67) F circuit = (g) ffW . 

The gates in G, labeled by the index [i, are described by unitary operators V^, so 
that the unitary operator describing the (noiseless) operation of the circuit is 

(68) Circuit = II % ' 

MSG 

where the product of operations is ordered in accordance with the definition of the 
circuit. The factors that appear in the multiplicand of l|68|) are in principle 
obtained from the fundamental solution to the appropriate underlying equation of 
motion, following the procedure outlined in Section 0] above. 10 

Each gate operates on a subset S M of the information elements, leaving the rest 
unaffected: 

(69) % = V^® (0/ H w) , 



It is extremely important to note that we are not discussing here the abstractly defined 
quantum computation itself, which is prescribed in advance, and is represented by the unitary 
operator U that appears explicitly in the second term under the norm symbol in the QCC given 
m G2J. Rather, we are discussing quantities that represent elements of a physical device that is 
to be used as an actual quantum computing machine. As such, the quantities under discussion 
here are to be regarded as "building blocks" for the first term under the norm symbol in 1121 . 
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where Ijji*) IS the identity operator on H^ l \ and V^ 1 * is the transformation effected 
by the nth gate, so that we may express the unitary operator describing the idealized 
circuit as 

(70) P-P = KhcuitPKircuit 

\neG j \neG j 

where the prime indicates that the second product reverses the order of the factors 
relative to the first product. 11 

Up to this point we have made no restrictions on the Hilbert spaces or the 
types of gates appearing in the specification of the circuit. We now describe three 
special cases of interest within the circuit-based paradigm: qubits, quantum con- 
tinuous variables, and liquid state NMR. We obtain the completely positive, trace 
preserving map P that characterizes the implementation of the circuit for each ex- 
ample, thus showing that the QCC provides the proper foundational presentation 
of a quantum computer for all the cases. 
(Case 1) Circuit-realization with qubits 

The great majority of the research in quantum computation has focussed on 
circuits for which the information elements arc qubits defined on a two dimensional 
Hilbert space C 2 , so that the full Hilbert space of the circuit is 

(72) ^circuit = (C 2 ) 8 ' ' . 

where \R\ is the cardinality of the set of computational qubits. 

Usually the set of gates is chosen from a relatively small set of operations, each 
of which acts on only 1 or 2 qubits. Specializing l|69|) to the case of a circuit built 
out of gates acting only on 1 or 2 qubits (this would be the proper description, for 
instance, of a machine that uses the universal set of quantum gates), we have 

(73) % = V^® I k , 
or 

(74) % = ® (g) h ■ 

where Ik is the identity operator acting on the copy of C 2 associated to the kth 
qubit. With these definitions, the operator P that characterizes the implementation 
of the qubit-based realization of a quantum computer designed according to the 
circuit-based paradigm is given by and JJDJ. We thus see that the circuit- 
based paradigm, on which a large amount of the research in the field is based, is 
properly described by the QCC. 

(Case 2) Circuit-realization with quantum continuous variables 

11 We strongly reiterate the message given in Footnote I1UI It is only coincidentally the case 
that the form of I7UI resembles that of Both the right- and left-hand sides of 1701 describe 
the physical device (i.e., P), not the abstractly-defined quantum computation U, and thus both 
sides of eq. l7UI arc to be associated with the first term under the norm symbol in the QCC given 
in 
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Alternatively, we may select different Hilbert spaces appropriate for quan- 
tum computation using quantum continuous variables (QCV). For instance, the 
ijW might describe the states of simple harmonic oscillators. The full Hilbert 
space of the circuit, and the gate operations out of which it is built, are then de- 
fined analogous to the above prescriptions for the qubit-based quantum computer. 
In this way we arrive at a specification of the quantum component P appropriate 
to the case of computation by QCV. Thus, the QCV realization of the circuit-based 
paradigm is also seen to be properly described by the QCC. 
(Case 3) Circuit-realization with NMR states 

The treatment of quantum computation by nuclear magnetic resonance (NMR) 
in liquids requires special consideration due to the fact that the NMR sample effec- 
tively contains many copies of the circuit carrying out the same computation. In 
this case we define the Hilbert space of the system as 



(75) Hnmr — H t 



circuit ' 



where N s is the number of copies of the circuit, which requires that we extend the 
definition of the unitary operation for the circuit as follows: 



(76) Vnmr = V£ 

We then obtain 



(77) Pp = VnmkPVumb. 

which describes the NMR-based quantum computer in the QCC. 12 

5.2.2. Quantum Circuits in the Presence of Decoherence and Dissipation. Thus 
far in this section we have restricted ourselves to a discussion of quantum circuits 
in the absence of decoherence and dissipation. This is reflected in our description 
of the gates as implementing purely unitary operations, as in 

JSHl - Even at this 

level of description the circuit is not necessarily error-free. Unitary errors derive 
from a situation in which the evolution of the quantum circuit is unitary, but the 
circuit does not implement exactly the desired unitary computation: 



(78) Vcu-cuit + U . 

Unitary errors can arise from either the design or the physical implementation of the 
circuit. Design errors arise, for example, due to the the fact that a universal set of 
quantum gates only allows for the implementation of an arbitrary unitary operation 
U to within an arbitrarily small tolerance |23j . In general there will then be some 
residual error implied by the very design of the circuit. Implementation errors 
result from inaccuracies in the physical parameters governing the unitary evolution 
associated with a gate as compared with the specification of those parameters by 
the design. For example, an interaction Hamiltonian may be applied for a longer 
time than specified, or there may be errors in the field strengths or couplings in the 
interaction Hamiltonian. 



'in connection with 1771 we repeat the admonition given in Footnote II II . 
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In addition to unitary errors, we also need to deal with errors resulting from 
decoherence and dissipation. We will refer to these simply as decoherent errors. 1 
In this case the gates are described by completely positive maps P^, where the 
hat on the P indicates that this quantum component is a single gate, and the 
index \i identifying the particular gate, as above. In addition to replacing unitary 
operations representing the gates by completely positive maps, it is also important 
to account for errors occurring in the transmission of quantum states from one gate 
to the next. This means that the transmission channels arc also represented by 
completely positive maps designated by the symbol P^. In effect, the transmission 
channels (including any quantum memories used to store the states) are regarded 
as gates that (ideally) implement the identity transformation: = I. If we call 
the set of transmission channels C, the completely positive map that describes the 
circuit is then 



(79) p = n h . 

iieGuc 

where the index \i now identifies both the transmission channels (in C) and the 
gates (in G). 

The operators P M are obtained by explicitly solving the equations of motion for 
the physical device that implements the gate. The result may be written formally 
as the sum of two terms, one of which represents the action of the unitary operator 
Vfj, describing the ideal gate as specified by the circuit design, and the other of which 
represents the effects of unitary implementation errors as well as decoherence and/or 
dissipation: 

(80) i> = (1 - e M ) (%ptf) + e^Q^p , 

where e M is the probability that an error occurs during the operation of the gate, and 
is a completely positive map that represents the effects of the error. Although 
we have indicated how this follows in principle from an explicit solution of the 
detailed dynamics of the gate (as described in jJJJ, an error model of this form 
is often assumed from the outset. The latter approach necessitates the choice of 
some specific operator to represent the errors. For instance, in investigating the 
properties of quantum error correcting codes one often invokes a "depolarization" 
qubit error model in which 

1 / 3 

(81) <3 M P= - p + ^VjPCTj 

V ^ 

where the aj are the Pauli matrices acting on a single qubit. The advantage of this 
approach is that it abstracts the physical implementation of the quantum computer 
from the design of the circuit while retaining the main features of decoherence that 
must be addressed in the development of any practical quantum computer. The 
disadvantage is that the abstraction must then be justified relative to the detailed 

^The case of liquid state NMR admits another type of error, which arises when the operators 
^ch-cuit appearing in 1761 effect different unitary errors on different copies (fc) of the circuit. These 
arc known as incoherent errors. 
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physical implementation of the computer in order to apply rigrously the results of 
any analysis based on (|H If) . 

Whether we arrive at l|8U|) by detailed calculation or by abstraction, we may 
use it in conjunction with (|T5|l to describe the operation of the entire circuit: 



(82) 



P-p 



Vr 



circuity V circuit 



nc 1 -^) 



e f Qfp , 



where Qf incorporates the effects of the individual gate errors. We define an error 
probability associated with the entire circuit: 



(83) C/ = l-JJ(l-e„) 

with which we have 



(84) P ■ p = (1 - e f ) V clrcuit pV c \ rcuit + e f Q fP . 

We will make use of this error model in applying the QCC to the problem of finding 
error thresholds for fault tolerance in <JH1 

We note that eqs. 1(80(1 and 1)84(1 connect theoretical analysis with experimental 
observations. A general theoretical method for obtaining an explicit expression 
for the quantum component, P, in terms of the underlying equation of motion 
for the system, has been given above in [Jl] From an experimental perspective, 
explicitly writing l(80|l or 1(84(1 is a goal of quantum process tomography, which is by 
definition the experimental method of determining the evolution of open quantum 
systems |33( . 

5.3. Adiabatic Quantum Computing Paradigm. We now show that the 
QCC also provides the proper framework in which to formulate the adiabatic quan- 
tum computing paradigm. In adiabatic quantum computing, the quantum compo- 
nent can be described by a Lindblad equation incorporating a Hamiltonian of the 
form: 



(85) ^adiabatic (t) = / (t) H + Q (t) H f 

where / and g are smooth functions of time with /(0) = 1, f(T) = 0, g(0) = and 
g(T) = 1, so that the adiabatic Hamiltonian goes smoothly from TCo to Ttf over 
time. The full Lindblad equation may then be written as 

where V(t) represents unitary errors and the terms involving Lj(t) represent inter- 
actions with the environment. 

The computation begins with an an initial preparation of the ground state of 
the Hamiltonian TIq. Provided the conditions for adiabatic evolution are satisfied, 
evolution under the exact adiabatic Hamiltonian takes the ground state of Tio to 



(86) —p(t) = -^[Wadiabatic(t)+V(i),p]+^ 
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the ground state of H / with a high degree of accuracy. We identify the operator U 
that appears in the QCC as a unitary operator that describes this desired behavior: 

(87) U : |0o) » \fh) ■ 

where |</>o) is the ground state of Ho, \ipo) is the ground state of Hf, and U is 
otherwise arbitrary. 

The physical realization of the adiabatic quantum computer is described by 
(|86[) . and thus implements the desired operation U only approximately. This is 
due to the approximation inherent in the adiabatic evolution itself, as well as any 
additional unitary errors, decoherence and dissipation. We express this fact quanti- 
tatively by using the fundamental solution of H86|l to derive the map P (as described 
above in jJJ that describes the operation of the adiabatic quantum component. The 
construction of the QCC then follows. 

Note that the special case of the adiabatic quantum computer is characterized 
under QCC(P,U, Ai enc , M.a ec , a) by two unique features: 

(1) the encoding and decoding maps, A4 cnc and Addec are identities, and 

(2) the QCC is required to hold only for the ground state of the initial Hamil- 
tonian, that is, for p = |</>o)(</>o| • 

5.4. Graph state-based (including cluster state-based) paradigm. We 

now consider the cluster-based approach to quantum computation f|25 | .|26 ) . \27\). 
which is formulated in terms of an array of two level quantum systems. The sys- 
tems in the array are referred to as sites and are elements of some set L which 
has a geometrical structure such as a 1 or 2 dimensional lattice; in addition to 
the geometrical structure there is also a flow structure which models the flow of 
information through the cluster. The two-dimensional Hilbert space corresponding 
to a site s £ L is denoted H s and the Hilbert space H of the entire cluster system 
is the tensor product 

(88) -^cluster = (g) Hs- 

seL 

A key concept in the cluster approach is measurement at a site s £ L considered as 
an operation on quantum mechanical states. As an operation on pure states, the 
measurement corresponds to a pair of self-adjoint projections E s , 1 — E s on H s . 
In terms of the cluster system, we identify the projection E s with a projection on 
-^duster defined by 

(89) E = E s ® (g) I Hm 

where Iu m is the identity operator acting on H m . The associated projective mea- 
surement on the cluster Hilbert space H c i uster is the completely positive map on 
density states 14 

(90) P E ■ p = EpE + (1 - E)p(l - E) . 

The cluster scheme is illustrated in Figure ^ 

A general class of cluster-type configurations associated to graphs has been in- 
troduced in the literature f|27). |1]). Given a graph Q — (nodes, edges), the lattice 
sites are the elements of nodes. Each lattice site a 6 nodes has associated with it 



We note that the message of Footnote II II applies to this equation. 
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Figure 1. Two-dimensional cluster configuration (the arrow de- 
notes time flow) 

a two-dimensional Hilbert space H a and an "entanglement" projection operator F 
that acts on the Hilbert space 

(91) ffgraph= (g H a 

aE nodes 

as follows: 15 To each (a, b) G edges, define the projector of the form 

(92) (g) lH n ®F latb) . 

nGnodes\{a,6} 

defined in terms of the Pauli matrices by 

(93) F M = 4 a) + 4 b) " ® 
Then define 

(94) F= [] ^(a.b)- 

(a,6)(Eedges 

f is a projection, since all the projectors F a o pairwise commute. 

5.4.1. The Cluster Measurement Scheme. In the cluster-based approach as ex- 
plained in |25j . a computation is realized by a sequence of projective measurements 
performed on different sites sei. The projective measurement that is to be per- 
formed at each step of the computation is determined by a scheme that specifies 
at which site to make the measurement and what observable to measure at that 
site. These two choices depend on the outcome of the preceding measurements. In 
order to specify this process, a discrete time "flow" between the sites is also given. 
This flow determines the sequencing of sites to measure. It is important to note 
however that the specific projective measurement taken at each site depends on the 
outcome of the measurement taken at the preceding site. 

At the level of specification, we can also consider the cluster measurement 
scheme as given by a multi-rooted tree T. In Figure |3 we illustrate such a tree 
which because of spatial limitations is singly rooted. The tree consists of nodes 
and directed branches. Each node on the tree T corresponds to a pair (s, A) where 
s G L is a cluster site and A is a two-level observable on H s . We will refer to s as 
the cluster site corresponding to the tree node (s, A). Each node I = (s, A) of the 



15 Such operators arc considered from the more general point of view as partial isometries 
below in Definition 15. 21 and in the discussion preceding it. 
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Figure 2. Example scheme of a cluster computation 



tree has two outgoing branches corresponding to the two possible outcomes of the 
measurement of A. These two branches correspond to the spectral projections of 
A, which wc denote E + (A), E~(A). 

It is important to note that many different nodes of T may correspond to the 
same cluster site, that is two distinct computation sequences may take us to the 
same node but at which different measurements are taken. In fact, in the usual 
cluster approach all computation sequences traverse the exact same cluster nodes. 
This means that at each horizontal level of the tree T, the nodes all have the same 
cluster site. In the general scheme outlined above, no such restriction exists. Thus 
we may consider schemes in which not only the subsequent measurement depends 
on previous outcomes, but in which the site at which the measurement is taken also 
dependent on previous outcomes. 

For clusters that are not arranged in a 1 dimensional array, we can avoid the use 
of multiply rooted trees if we allow cluster systems that are not restricted to two- 
level systems. This means that the Hilbert space H s corresponding to a site seL 
can have arbitrary dimension. In that case, the corresponding node observables are 
allowed to have more than two outcomes and in particular, the specification tree 
may not be binary. 

A complete path (that is one which originates at the root node and terminates 
at a leaf node) of T is specified by a sequence of measurement outcomes, where 
the observable measured is associated to a node of the tree. Mathematically, each 
measurement outcome is expressed by one of the spectral projections associated to 
the measurement and graphically represented by an edge of the tree. A complete 
measurement outcome associated to a path is a sequence of projections, where each 
projection corresponds to a traversal of an edge of the tree as follows: 



(95) 





The projective measurement associated to the cluster consists of the projective 
measurements on sites following the branches of the cluster scheme tree. We can 
explicitly write this down: 
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Theorem 5.1. The projective measurement on a cluster is of the form 

(96) P T - P = (li E -) p(li E -) 

PePathTVe'P ' Vep ' 

where E T is a projector on the cluster site corresponding to the edge of the cluster 
scheme tree. Note that for each V G Path(T), all the E T with r G V commute. 

5.4.2. Quantum Components in the graph state-based paradigm. We now show 
that the graph state-based paradigm of quantum computation is properly described 
by the QCC. The paradigm has been described in the literature as employing an 
entangled substrate upon which a series of conditional projective measurements is 
performed. The projective measurements are used to carry out a computational 
algorithm, but can also be used to read input from a macroscopically observable 
input register or output to a macroscopically observable output register. The en- 
tangled substrate corresponds to some particular Hilbert subspace of entangled 
vectors which may be characterized in various ways, for instance as the range of an 
entangling operation or as the solutions to some eigenvalue equations. 

Correspondingly, we should expect that the mathematical formulation of the 
graph model of a quantum component also consist of two parts (We assume as given 
the Hilbert space i? gra ph): 

(1) An initial "entanglement producing" operation of some kind. 

(2) The projective measurement on i? gra ph corresponding to the graph scheme. 

In our approach, we have already discussed how to formalize the step (2) above. 
However, there are several choices for the entanglement generating step (1). In the 
examples discussed in the literature, these are given by a self-adjoint projection 
operator, but it is natural from our viewpoint to consider more generally partial 
isometries V : H — ► H . 

Definition 5.2. A graph state-based quantum component C is a pair (V, T) 
given by a cluster scheme T and an "entanglement" partial isometry V on i? gra ph- 
The completely positive map associated to C is defined by 

(97) Qe-p = E V V P V^E V 

•PGPath(T) 

Thus, we see that, just as for the circuit-based paradigm and the adiabatic 
quantum computing paradigm, the QCC provides an over-arching framework in 
which to formulate the graph-state based paradigm of quantum computing. More- 
over, in this section we have generalized the definition of graph state- (and cluster 
state-) based quantum computers that has previously appeared in the literature. 
Our generalization consists of two features: (1) our definition allows for arbitrarily 
different measurements to be carried out at different nodes at the same level of 
the multi-rooted tree T, and (2) our definition replaces the use of a self-adjoint 
projection operator as an entanglement generator with the more general notion of 
a partial isometry (which includes self-adjoint projections as a special case). 

6. Error Thresholds and Fault Tolerance 

6.1. Introduction. In i i2.3.3l wc showed that the recently discovered approach 
to error correction known as "operator quantum error correction" is in fact a special 
case of the general QCC formulation. Given the general applicability of the QCC to 



A THEORY OF PHYSICAL QUANTUM COMPUTATION 



31 



all quantum computing paradigms (cf flBjl. as well as to all techniques for protection 
against errors (including quantum error correction, decoherence-free subspaces and 
noiseless subsystems), the QCC thus provides a unified framework for a fully general 
analysis of fault tolerance in quantum computing. We refer to this as operator 
quantum fault tolerance (OQFT). 

As an example of OQFT, in this section we describe the application of the QCC 
to the analysis of error thresholds for fault tolerance in the circuit paradigm. To 
make contact with previous research, we begin by discussing this subject from the 
perspective of the well established method based on the analysis of error probabil- 
ities [28], El; EH, EU. HH, HI- Since the Q cc in fact provides the underlying 
framework in which to study any issues associated with physical quantum compu- 
tation, we then reformulate the problem in terms of the QCC. This allows us to 
relate the results of the two approaches, and to determine the extent to which the 
previously utilized approach (based on the method of error probabilities) is in fact 
justified based on the insight provided with the QCC. 

6.2. The Method of Error Probabilities. In this section we briefly de- 
scribe the method of error probabilities. We begin by identifying a quantum oper- 
ation that we wish to implement and specifying a circuit that (ideally) implements 
the operation. We then identify an error model for the gates in the circuit that ac- 
counts for the inevitable effects of dissipation and decoherence that come into play 
when the gates are implemented with real devices. By hypothesis, the probability 
that an error occurs in this "direct" implementation of the operation is too high 
for it to be useful as a component in a quantum computer. 

In order to render the circuit more fault tolerant, we specify a second, more 
complicated, circuit using quantum error correction. The circuit now operates on 
the set of encoded qubits obtained by encoding the logical qubits of the direct 
implementation using a QECC. The gates in the original circuit are replaced by 
collections of gates that operate on the encoded qubits. Additional gates are added 
to carry out the QECC's recovery procedure wherever an error is detected. The 
error model is then applied to the gates in this new circuit. The error correction 
clearly provides some degree of protection against errors, but it also necessitates a 
larger number of gates, so that there are then more opportunities for errors to occur. 
In order to determine whether this procedure has improved the fault tolerance, we 
compare the probability of an uncorrected error occurring during operation of the 
QECC version (e^ E ) with the probability of any error occurring in the "direct" 
implementation (e^ lrect ). If the QECC error probability is smaller, that is, if 

QECC 

( 98 ) iiSF < 1 - 

then the procedure has been at least partially successful. If the error probability 
of the new circuit is still too high, we can reduce the error probability further by 
concatenating the code, that is by encoding the qubits used in the QECC version 
using the same QECC and by redesigning the circuit to handle the second level of 
encoded qubits. By repeating the process of concatenation we can arrange that the 
error probability be made arbitrarily small. 

The key point is that 1)98(1 can be shown to hold only if the failure probabil- 
ities of the individual gates are less than some threshold value that depends on 
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the relative complexity of the encoded vs. un-encoded versions of the circuit. The 
thresholds appearing in these conditions are known as "error thresholds." If the 
threshold conditions are satisfied, then the use of concatenated codes will provide 
fault tolerant operation to within some specified tolerance. The problem of achiev- 
ing fault tolerant quantum computation is reduced to the problem of constructing 
implementations of the primitive gates that satisfy the error threshold conditions. 

6.3. Operator Quantum Fault Tolerance and the QCC. We now recon- 
sider the above problem from the perspective of OQFT by making use of the QCC. 
Our goal is to identify a quantum component P, that implements (approximately) 
a quantum computation, U, that is, we write the quantum computer condition, 
QCC{P,U,M cnc ,M dec i C^J; 

(99) \\M dcc {P ■ (M cnc (p))) ~ UpU% < a 

for some suitable choice of encoding and decoding maps, M cnc and Md cc . For 
simplicity of notation we define an operator 

(100) P = M dcc ■ P ■ M cnc , 
so that the QCC becomes 



(101) \\Pp-UpU*\\i < a . 

We begin by developing a "zero-th order" implementation that does not use 
a QECC. In the absence of errors, we assume that the implementation faithfully 
implements the computation U: 



(102) 



p(0) 



P = M 



dec 



V- 



(0) 



(AW) V c 



(o)t 

circuit 



UpUl 



With the error model l|84|l we have 



(103) 

\\P (a) p-U P U^\\ = 

where, for generality, we have subsumed the encoding and decoding maps into the 
definition of Qf. (The maps arc identities for the zero-th order circuit.) With i|lU2|) 
the left hand side of the QCC takes the simple form 



1 - e ( f 0) ) M dec 



circuit (M e 



(0)t 
cuit 



e^Q^p-Uptf 



(104) \\pWp-Uptf\\ =ef\\Qfp-U P U^\\ l . 

Note from l|83|l that the failure probability for this circuit is, to lowest order, linear 
in the error probabilities for the gates which make up the circuit: 



(105) 4 0) ~om, 

(recall that e M represents gate error, cf (HSJ) a fact which, as we shall see, is crucial 
to the derivation of an error threshold. 
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We next construct the "first order" implementation of U, which operates in 
the codespace of the QECC. By the preceding arguments, this implementation will 
satisfy 



(106) \\P^p-UpUm =ef\\Qfp-UpU% . 

Following |24j . we consider the case that errors affect the qubits in the circuit 
independently and that the QECC recovery procedure is sufficient to correct a 
single error in any one qubit. In that case, the encoded circuit will exhibit an 
unrecoverable error only if two or more single qubit errors occur. In this case, the 
error probability will be quadratic in e M to leading order: 



(107) ef-OM 2 . 

At this point, the QECC has not eliminated all errors, but has resulted in a cir- 
cuit with error probabilities that are quadratic, rather than linear, in the error 
probabilities of the primitive operations. 

We have now introduced two implementations of the computation. If cither 
implementation satisfies the QCC, then there is no need to continue. If the imple- 
mentations do not satisfy the QCC, then it is meaningful to ask whether this can 
be achieved by concatenating the code. In order to answer this question, we begin 
by asking another, related question: has the QECC improved the fault tolerance of 
the implementation relative to the QCC? In other words, we wish to know under 
what conditions it is true that for all p 



(108) ||P«p - UpU% < ||P(°) P - Uprf\\i 

or alternatively 16 



\\m P -u P u% 

||p(o; 

Using HI04jl and this becomes 



(109) sup„ " , - — — < 1 

P \\PWp-UpU% 



4 X) \\Q?p-UpU% 1 

(110) sup pl5r^fe — ^77tT <1 - 

e} WQyp-UpWW! 

The above equation represents the extension of to OQFT obtained by using the 
general framework provided by the QCC. Clearly it does not reduce to a simple ratio 
of error probabilities, as in (|98|) above. The reason for this is that this formulation 
of the error threshold problem based on the QCC takes into account not only 
error probabilities, e/, but also expressions involving the norms of operators that 
characterize the "strength" of the errors. To see this, note that we began by 
assuming an error model represented by the operation for the zcro-th order 
implementation. The error model in the encoded implementation is represented, in 
general, by a different operation, There is no reason to suppose that these 



16 Rclation I1U9I is actually a stronger condition than I1U8I for infinite-dimensional vector 
spaces. 



:S4 



GERALD GILBERT, MICHAEL HAMRICK AND F. JAVIER THAYER 



error models are equally effective in perturbing the computation. This is reflected 
in the difference in the norms: 



(111) \\Qfp-UpU%^\\Qfp-UpU% . 

To further emphasize this point, we obtain the above result l|98|) based on the 
assumption that the error models are "commensurate" in the sense that 

(112) \\Qfp-UpU%*\\Qfp-UpU% • 

We shall shortly return to the question of whether this is a good approximation. 
With this assumption, the condition l|110|) then becomes 

e (1) 

which is identical to the result (|98|l obtained by the method of error probabilities. 
We obtain the form of the error threshold by noting from l|105|) and (|107|l that 
the numerator and denominator of H113(l are, to lowest order, quadratic and linear, 
respectively, in e^, so that 

V A 



(114) 



At this point, it is straightforward to obtain a threshold if we set the e M equal to 
each other, and take e = e^. Then one obtains the error threshold comparable to 
those obtained in , [13] , pS] , as 

(H5) £ <JP^. 

The desired behavior of the concatenated QECC described above follows if we 
successively apply the approximation (|112|1 at each level of concatenation: 



(116) \\Qf x \-XJptf\\ X K\\Qfp- UpU% . 

We note that Aliferis, Gottesman and Preskill |2j also relate the ratio of error 
probabilities for successive levels of concatenation to an overall measure of the 
"accuracy" of the quantum computation. Their approach differs from the one 
described here in three important ways: 

(Contrast 1) Accuracy in [5] is defined in terms of the probabilities pi of the out- 
comes i of measurements on the final output state for an ideal, as compared with 
a noisy, circuit: 

(117) £br isy -M dcal i • 

i 

In contrast, we define the accuracy of the implementation by the QCC. 
(Contrast 2) The threshold proofs in |2] rely on proofs that the implementations 
at each level of concatenation are conditionally correct relative to the noise model. 
In this case the ratio of error probabilities e( l+1 )/eW is automatically the quantity 
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of interest in comparing performance at each level of concatenation. Here, in con- 
trast, the QCC provides the figure of merit at each level of concatenation, and the 
dependence on error probabilities is inferred. 

(Contrast 3) As a consequence of the two preceding points, |2j makes contact with 
the notion of accuracy only at the highest level of concatenation, at which the entire 
quantum component may be viewed as a black box. Here, the QCC is applied 
systematically at each level of concatenation. 

At this point we have shown that we can obtain the standard form of the error 
threshold result from the QCC by introducing the assumption (|116|l on the relative 
strengths of the error models appropriate to successive levels of concatenation of the 
QECC. We now investigate the validity of the approximation. We begin by making 
some reasonable assumptions about the error operators and the initial state of the 
computer. We note that the operators are trace preserving, and thus describe 
the evolution of the quantum component if a failure has in fact occurred, as can 
be seen by setting ej = 1 in l|84|) . It is then reasonable to expect that the state 
resulting from its operation on p will be "close" to the maximum entropy state 17 
in the sense that, for small 5 Pl 



(118) 



PI 



\ P < 5 P 



where pi is the maximum entropy state, and the value of p identifies the Schatten 
p-norm associated to the corresponding Schatten p-class (cf l|126|l ). 
On the other hand, it is normally the case that the input state for the quantum 
computation is a pure state, and thus so is the state UpW. In this case, it is 
straightforward to show that 



(119) 
and 



\ Pl -UpU%=2-- 



(120) \\pi-UprfWac 

where N is the dimension of i?i og icai- Since 



1 



(121) \\Qfp-UpUH = || \Qfp-pi] + [pi-Uprf] || , 



we have, by triangle inequalities, 
2 



(122) 
or 



-^-S 1 <\\Qfp-UpU^\\ 1 <2 



N 



Si 



(123) 



1 

TV 



< 



SCO 



l^Note that this is not a good assumption for errors that operate locally on only one qubit 
in a larger set of qubits, leaving the others unaffected. This important case is a topic for further 
study. 
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where we have assumed implicitly that N 3> 1 and S p <C 1 . Since these expressions 
hold for any value of i, (|116J) is a reasonable approximation with either choice of 
norm under the conditions that 

(1) (|118|l holds for some (5 p < 1, and 

(2) the initial state of the computation is a pure state. 

7. Conclusions 

In this paper we have presented a fundamental, unifying framework for describ- 
ing physically-realizable quantum computing machines. This is concisely stated in 
the form of the Quantum Computer Condition (QCC), an inequality that incor- 
porates a complete specification of the full dissipative, decohering dynamics of the 
actual, practical device used as the quantum computing machine, a specification of 
the ideally-defined quantum computation intended to be performed by the machine, 
and a quantitative criterion for the accuracy with which the computation must be 
executed. 

With the QCC we prove the fundamental Encoding No-Go Theorem that iden- 
tifies the amount of damping (including dissipative and decohering effects) for which 
physically-realizable fault-tolerant quantum computing is not possible. We provide 
a rigorous definition of damping, and explicitly calculate a universal critical damp- 
ing value for fault-tolerant quantum computation. This theorem can be used in 
principle to solve practical problems involving quantum computer design. 

In this paper we have also presented an existence proof for fundamental solu- 
tions to useful classes of time- dependent generalizations of the Lindblad equation. 
This can provide a useful tool in analyzing a wide variety of open quantum me- 
chanical systems. 

We have demonstrated that the entire formalism of operator quantum error 
correction (OQEC) can be obtained from the QCC as a special case. By allowing 
for the possibility of residual errors, the general formalism of the QCC enables us 
to generalize OQEC to "operator quantum fault tolerance" (OQFT). Since we have 
demonstrated that OQEC is in fact a particular reduction of the QCC, and since 
standard quantum error correction (QECC), decoherence-free subspaces (DFS) and 
noiseless subsystems are all special cases of OQEC, we have discovered that QCC 
applies in general across all these approaches. 

As an initial application of the OQFT concept, we have begun the exploration of 
the application of QCC to the problem of establishing thresholds for fault-tolerant 
quantum computation by showing that the standard approaches to this problem 
can be motivated within the framework of the QCC. 

Research in quantum information science has resulted in the discovery of seem- 
ingly different paradigms for quantum computation, including the circuit-based par- 
adigm, graph state-based paradigm and adiabatic quantum computing paradigm. 
In this paper we have explicitly demonstrated that these paradigms are not in fact 
distinct at a fundamental level, but arc all dcscribablc within the unifying frame- 
work provided by the QCC. In the particular case of the graph state-based paradigm 
(which includes cluster state-based models), we not only show that the paradigm is 
a manifestation of the unifying picture provided by the QCC, but also introduce a 
definition of graph state-based quantum computers that generalizes the graph state 
models previously defined in the literature. 
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Future work motivated by our results should include applications of the Encod- 
ing No-Go Theorem to diverse problems pertaining to practical quantum computer 
design and implementation. It would also be of interest to further explore the 
physics of the operator quantum fault tolerance (OQFT) generalization of OQEC 
presented in this paper. Specific work along theses lines should include further 
application of the QCC to obtaining error thresholds for fault-tolerant implemen- 
tations of quantum computers. It would also be fruitful to explore applications of 
the quantum computer condition to situations in which quantum process tomog- 
raphy techniques are used to experimentally characterize the quantum component 
described by the completely positive map that appears in the QCC. 
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Appendix A. Banach Spaces of Operators 

A linear map T on a Banach space is a contraction iff its norm is < 1 . 

Let H be a separable Hilbcrt space. L(iJ) denotes the space of bounded oper- 
ators on H with the operator norm, K(i?) denotes the normed closed subspace of 
compact operators of li(H). In case H is finite dimensional, these spaces arc iden- 
tical. Assume now H is infinite dimensional; we consider other Banach spaces of 
compact operators defined by eigenvalue decay conditions and whose norms reflect 
the rate of decay of the eigenvalues. Specifically, let T S K(7J), then \T\ = vT^T 
is a non-negative compact operator and so by the spectral theorem, has a complete 
set of eigenvectors with eigenvalues that can be ordered in a sequence 

(124) so(T)> S i(T)>--->.s n (T)>0 

which converges to 0. The following properties are well-known (see [12| : also [7] on 
which this discussion is based). 

T(H) is the Banach space of trace-class operators T on H with the norm 

oo 

(125) imii = x>coi 

k=0 

More generally, the Schatten p-class T p (H) is defined by the condition 

{OO i 
5>(iy 
fc=0 ' 

with the norm given by || • The operator norm for any T £ L(_ff), denoted by 
\\TWoo, is defined as the supremum of ||Ta;|| for x e H of norm < 1. If T £ K(H), 
the operator norm is also the supremum of the eigenvalues of \T\. 

Appendix B. Completely Positive Maps and Fundamental Solutions 

In the following, H denotes a complex Hilbcrt space. We will consider various 
Banach spaces of bounded operators on H : these are discussed in Appendix lAl 
above. 

We will consider completely positive maps on both ~L{H) on the Schattcn classes 
T P (H) and especially on the trace-class operators T(H) = Ti(H). 
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The basic fact about completely positive maps we use is the Kraus representa- 
tion. 



(127) P(p) = '£x i pX} 



Proposition B.l. A P : T(H) — ► T(H) is a completely positive contraction 
if and only if it is of the form 

iei 

where Xi £ L(iJ) with 

(128) Y. X i Xi ^ 1 

i 

We will consider only trace preserving completely positive maps P, that is which 
satisfy 

(129) tr(PQo)) = tr(p) 

Proposition B.2. Suppose P is a completely positive map given by the Krauss 
form (tT27jl . 

(1) A necessary and sufficient condition P be trace-preserving is that 

(130) E^' = l- 

iei 

(2) A necessary and sufficient condition that P be bounded in the operator 
norm is that 

(131) ^XiXjeLiH). 

iei 

Note that if T is a completely positive trace-preserving and operator norm 
continuous map, then by interpolation T is also norm continuous on the Schatten 
p-classes. 

Proposition B.3. If P is a completely positive trace-preserving map, then the 
adjoint of P on Jj(H) defined by 

(132) tr(P t (T)p)=ti(TP(p)) 

is a unit preserving completely positive map L(if) — > ~L{H). Its Kraus representa- 
tion is 

(133) P\T) = Y^,X\TXi 

iei 

Note that the completely positive unit preserving maps L(i?) — » L(iJ) which 
are adjoints of completely positive trace-preserving maps on T(H) — > T(iZ) can 
be characterized precisely as those which are continuous mappings L(if) — > ~L{H), 
where L(iJ) has the ultraweak topology. 

B.l. Completely Positive Semigroups on T(H). We need to first estab- 
lish that each generalized Lindblad-type operator A(t) (cf eas. l|S5|l and (HHU) gen- 
erates a semigroup of completely positive contractions on T(H) with respect to the 
trace- class norm \\ ■ ||i. Wc will also consider boundcdness properties relative to the 
operator norm || ■ Hoc, although in general the corresponding semigroups may not 
be contraction semigroups in this norm. 
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We begin with a general result characterizing infinitesimal generators of strongly 
continuous positive (or completely positive) semigroups on the Banach space T(H). 
Recall that a one-parameter semigroup {T t } t >o on a Banach space E is said to be 
of class Co iff for every x £ E, lim t _» s T t {x) = T s {x). If {T t } t >o is a Co-semigroup, 
then there arc positive constants M and (3 such that 

(134) ||T t || < Me tp . 
Moreover, 

(135) Ax = lim /i _1 (T ft a; - x) 

is a densely-defined operator called the infinitesimal generator of {T t } t >o 

Theorem B.4. Suppose A is a densely-defined operator on T(H) which gener- 
ates a contractive semigroup {T t }t>o on T(_ff). A necessary and sufficient condition 
the operators {T t } t> o be positive (respectively completely positive) is that for each 
A > 

(136) R(A, A) = (XI - A)' 1 

( which is defined by the Hille- Yosida Theorem ) be positive ( respectively completely 
positive). The operators T t are trace-preserving iff in addition for all A > 0, 

(137) tr(R(A, A)p) = A -1 tr(p) 

for every p G T(H). 

The family of operators {T t }t>o extends to a Co-semigroup on L(if) iff there 
are constants M' and (3 1 such that for all A > /3' and for all p G T(i7) and non- 
negative integers m: 

(138) || R(A,^4) m /9||oo < M'(A- /3')~ m |l/°iloc- 
In this case, we have the explicit bound 

(139) \\TtWoo < M'e t(3 ' Vt>0. 

Remark B.5. It suffices that the property (|137J) hold for density states p. 

Proof. To avoid duplication, we refer only to the assertions for complete pos- 
itivity. By the general Hille- Yosida theory, if A is an infinitesimal generator of a 
contractive semigroup on T(H), the resolvents 

(140) R(\,A) = (A- A)" 1 
are defined for all A > and by Theorem 3.1.3 of |30j . 

poo 

(141) R(A, A) = / e" xt T t dt. 

Jo 

i.e., the resolvent is the Laplace transform of {T t }t>o- In particular, if T t consist 
of completely positive operators, then the resolvent operators are all completely 
positive. 

Conversely, let 

(142) A\ = A(AR(A, A) - I) 
Then for each t > 0, 

(143) exp(L4 A ) = e.- xt cxp (t\ 2 R(A, A)) 
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which is clearly completely positive and it is known that for each t > 0, 

(144) T t = lim cxp(tA x ) 

X — >oc 

in the strong operator topology. Thus T t is completely positive. 

To deal with the trace preservation properties of T t , note that if S is a bounded 
operator on T(iJ) for which 

(145) tv(Sp) = atv(p) 
then 

(146) tr(e s p) =E tr ^ =Eff tr ^) = e ° tr ^ 

fc=0 ^ ' ' k=0 

Thus, 

(147) tr(expiA Aj o) = e~ xt e t>? 1/A tr(p) = tr(p). 

By <|144|> . it follows that T t is also trace preserving. Conversely, if T t is trace 
preserving, 

/>oo />oo 

(148) tr(R(A, A)p) = \ e~ xt tr (T t p)dt = / e _At tr(p)dt = A -1 tr(p). 

Jo Jo 

□ 

B.l.l. Examples of Completely Positive Semigroups. Our analysis of the gener- 
alized Lindblad equation will reduce to an analysis of the solution in two important 
cases: 

(Case 1) Unitary Evolution. In particular, if TL is a self-adjoint operator on a 
Hilbcrt space H , then the family of completely positive mappings 

(149) P«(p) =e im pe- im 

is a one-parameter group of completely positive maps. Its generator on the trace- 
class operators is formally given by the operator 

(150) P»i[H,p]. 

This expression is only formal, because it is not defined for all p. Nevertheless, the 
infinitesimal generator is densely defined on the space of trace-class operators and 
it is an extension of (|150|) for the finite rank operators on the domain of TL. 

(Case 2) Dissipative Operators. Another type of infinitesimal generator we will 
consider are operators of the form 

(151) Cp = J2\L j pL]-±{L]L j ,p} 

3 

where braces denote the anti-commutator. 
Lemma B.6. Suppose 

(152) Y^l-,l-, " •'• 
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Then the operator given by l|151|l is bounded on T(iJ). If in addition 

(153) »,/..; L;//i - 

3 

then C is a bounded operator on ~L(H). 

PROOF. Let C be the operator norm of ^jL^Lj. To show the map C is 
defined and continuous on T(H), it suffices to show Cq : p ^ J^j LjpLij is defined 
and continuous on T(H). However, if p > 0, 

(154) tr(£ (p))=^tr(L,pi}) 

i 

(155) =X>(p4^) 

o 

(156) =^{p 1/2 L)L jP ^ 

(157) = tr(V /2 (£4^V /2 ) < Ctr(p) = C||p|| x 



thus for arbitrary self-adjoint p, 

(158) \\£o(p)\\ 1 = \\C (p+-p-)\\ 1 

(159) < ||£ (p + )||i + ||/:o(^)||i 



(160) <c(j|p + ||i + ||p-||ij =cy| 

If UnSJ, suppose < T < 1: 

(161) < l ] tl j < L ) L 3 < C1 H 

j j 

Thus, for arbitrary T, 

(162) iirTii^ii^L^jii^ + qiTiu 



□ 



It was established by Lindblad |21| (and not too hard to show directly) that 
if l|152|) holds, the £ generates a uniformly continuous a completely positive semi- 
group (relative to the operator norm on T(H)). In this case the semigroup is given 

by 

(163) et£ = E^- 

fc=0 

B.2. Perturbation of Completely Positive Generators. In order to show 
that the generalized Lindblad operators given in equation l|60|l arc completely pos- 
itive generators, we need to establish a perturbation result analogous to the Kato- 
Rellich theorem. 

A linear map B on T(7?) is trace annihilating iff 

(164) tr(Bp) = 
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for all p € T(H). For example, an operator of the form (|151|l is easily seen to be 
trace annihilating. 

Using the Trotter-Kato product formula [5]), we can show that generators 

of completely positive contractive semigroups have a sum which is also a generator 
of a contractive semigroup, provided the sum generates a contractive semigroup. 

Proposition B.7. Suppose B is a bounded operator of the form (|151(l . If A 

is a generator of a completely positive semigroup then so is A + B . . 

COROLLARY B.8. The generalized Lindblad operators given in eguation l|6U|) 
generate a completelty positive semigroup of contractions on T(H). 

We now extend the results of Lindblad and Davies to allow for time varying 
Hamiltonians by relying on results of Kato. The need for this arises since in some 
circuit-based models the various gates arc implemented by varying the Hamiltonian 
(see for instance [23j §7.7.2). We make the assumption that the dissipative effects 
are bounded which simplifies the analysis considerably. 

B.3. Solving the Generalized Lindblad Equation. In some cases it is 
possible to solve the generalized Lindblad equation [22) . However, by solution 
we mean an expression for the fundamental solution P t . s as a limit of product 
of exponentials. Though this expression will almost never provide a closed form 
solution, it will provide enough information to obtain an estimate of how well a 
unitary (or partial isometry) can be implemented by one of the operators Pt. s - The 
two tools we use are the Trotter-Kato product formula and the explicit form of the 
solution of a time-dependent equation as a time ordered product of exponentials 
given in the proof of §4.2 of |30| . 

A precise formulation of a set of conditions which guarantees the convergence 
of the products in the next two theorems is given in Theorem IB. Ill These results 
comprised by Theorems IB. 91 and IB. 101 are restatements of assertions contained in 
the proofs in §4.2 of [371]. 

Theorem B.9. Under suitable conditions, the fundamental solution P t:S for (|62|) 
is given by 

n-i 

(165) P t .s = str-lim A ^o J| cxp((r fe+ i - r k )A(r k )), 

where s = Tq < r\ < • • • < r„_i < t and max \r k +\ — r k \ < A. Each P t:S is 
completely positive and trace preserving. 

Theorem B.10. Under the same assumptions as the previous theorem \B.fA 

(166) exp sA(t) = str-lim„^oo cxp ( —C(t) ) exp ( —Hit) ) 

\n J \n J 

B.4. Existence of Solutions. We will restrict our attention to bounded time 
varying perturbations of a fixed self-adjoint operator acting on T(H) via a com- 
mutator as in l|167|) below. The result we state is not the most general possible, 
and the early results of Kato [14| suffice for its proof. We follow the treatment in 
Chapter XIV, §4 of |34) which is a more readily available reference. 

Theorem B.ll. Suppose Ji is a self-adjoint operator, {B(t)} tG [ oo [, {Lj(t)} te [ oo [, 
1 < j < n are families of bounded operators, all of which are continuously norm 
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differentiable as a functions oft, then there is a fundamental solution Pt, s for Ij62(l 
where 

n , i 
(167) A(t)p = -i[H + B(t),p] +J2 [ L i(t)pL](t) - -{L){t)L 3 (i) lP } 



j'=i x 

The solution is a given by a limit of a time-ordered product of exponentials <|165|) . 

Proof. There are various technical assumptions for a family A{t) of operators 
that need to be checked in order to apply Kato's Theorem. The first of these is the 
independence of dovn A{t) of the parameter t. Under our assumptions 

(168) A(t)=A + C(t) 

where C(t) : T(H) — > T(H) are bounded operators and A is the infinitcsinal 
generator of a contractive semigroup. Indeed, 

(169) Ap = -i[H,p] 

is the infinitesimal generator of a group on T(7J) and 

(170) C(t)p = -i[B{t),p] + J2 {Li{t)pL){t) - \{L\{t)Li{t), p}^ . 

is by assumption a bounded operator on T(H). In particular, all the operators A(t) 
have the same domain dom(A). 

We now address the remaining assumptions in Kato's theorem. For any A > 0, 



(171) A - A(t) =X — A- C(t) = (I + C(t) R(A, A))(A - A) 

For A sufficiently large R(A, A)C(t) has norm < 1, so the Neumann (geometric) 
series for inverses (see [5], Chapter VIII, §3) I + R(A, A)C{t) is invertible. Thus we 
can write, 

(172) (A - A^))- 1 = R(A, A)(I + R(A, A)C(t))" 1 
Thus, 

(173) B(t, s) = (A - A(t))(\ - A(s))- 1 = (I + C(t) R(A, A))(I + C(s) R(A, A))- 1 

is well defined and by our assumptions B(t, s) is a norm differentiable function 
jointly in the variables t, s. This implies the remaining conditions in the hypothesis 
of Kato's theorem. It only remains to observe that presence of the parameter A, 
instead of 1 as actually stated in Kato's theorem is immaterial, since solutions of 
equations 

(174) j t p(t)=A(t)p(t) 



are trivially affected by adding a constant scalar to A(t). 



□ 
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Appendix C. Solving the Ersatz Quantum Computer Condition 

Consider the ersatz quantum computer condition (£QCC) given in 

(175) P-p = UpU f . 

Note that this can be obtained by setting the encoding and decoding maps to 
unity, and setting a = in the QCC given in l|12|) . A simple observation shows 
that "solving" (|175f) is completely equivalent to obtaining the noise-free part of a 
communication channel. 

Theorem C.l. Suppose H is finite- dimensional. If P is given by the Kraus 
representation (|127l) . given a unitary U , the set of p G L(i?) satisfying (1175(1 is the 
*-subalgebra o/L(iJ) given by 

(176) % P>U = {pe L(H) : V* G J, [p, U ] X,] = 0} 

Proof. The solutions of (|175fl are the fixed points of the completely positive 
map Q defined by the equation 

(177) Q. p = Y J U^XipXtu. 

iei 

Now apply jEj, Theorem 2.1. □ 

In the above theorem, the assumption H is finite-dimensional is essential, 
see P2- Since H is finite dimensional 21 is an algebraic direct sum of algebras 
isomorphic to full matrix algebras. 

Proposition C.2. Let {E k } k€ j be the set of finite- dimensional minimal central 
projections of 21. Then 

(178) 2t K = E K f&E K 

is an algebra of operators on the range H K of E K ( which is a finite dimensional 
space). It is isomorphic to a full matrix- algebra of finite multiplicity. 
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